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' Abstract 

in 

We use techniques from soft-collinear effective theory (SCET) to derive renormaliza- 
tion-group improved predictions for single-particle inclusive (1PI) observables in top- 
quark pair production at hadron colliders. In particular, we study the top-quark trans- 
verse-momentum and rapidity distributions, the forward-backward asymmetry at the 
Tevatron, and the total cross section at NLO+NNLL order in resummed perturbation 
theory and at approximate NNLO in fixed order. We also perform a detailed analysis 
of power corrections to the leading terms in the threshold expansion of the partonic 
hard-scattering kernels. We conclude that, although the threshold expansion in 1PI 
kinematics is susceptible to numerically significant power corrections, its predictions for 
the total cross section are in good agreement with those obtained by integrating the top- 
pair invariant-mass distribution in pair invariant-mass kinematics, as long as a certain 
set of sub leading terms appearing naturally within the SCET formalism is included. 



1 Introduction 



The top quark is the heaviest elementary particle known to date. Since its mass m t = (173.3 ± 
1.1) GeV [1] is on the order of the electroweak scale, the top quark is a crucial tool in the study 
of the electroweak symmetry breaking mechanism. The top-quark mass is an important input 
parameter in electroweak fits [2], and plays a role in the investigation of many Beyond the 
Standard Model (BSM) scenarios. The dominant top-quark production mechanism at hadron 
colliders is the simultaneous production of a top-anti-top pair. Several differential distributions 
related to the production of top-quark pairs, such as the tt invariant-mass distribution and 
the transverse-momentum distributions for the top quark, allow to investigate the existence of 
possible heavy s-channel resonances [3-6] predicted by many BSM scenarios. Moreover, the 
rapidity distribution of the top quark can be used to directly calculate the forward-backward 
asymmetry at the Tevatron [7-9], an observable of much interest because of its potential 
sensitivity to new physics. 

So far, the best measurements of the top-quark mass, couplings, and production cross 
sections have been performed using Tevatron data, where this particle was discovered in 1995 
[10]. In the last 15 years a few thousand top-quark pair events were studied by the CDF and 
DO collaborations. The LHC is now running at a center-of-mass energy of 7 TeV, and very 
recently the first measurements of the total pair-production cross section was performed by 
the CMS [11] and ATLAS [12] collaborations. Hopes are that it will be possible to obtain an 
integrated luminosity of up to lfb" 1 , which would produce roughly 150,000 tt events before 
selection [13]. With the planned increases in the center-of-mass energy and luminosity at the 
LHC, it will eventually be possible to observe millions of top quarks per year. The increase 
of the production rate will induce a decrease of the experimental errors on the measured 
observables. For the total inclusive cross section the relative experimental uncertainty is 
expected to become of the order of 5%-10% [14]. 

In view of the precise measurements at the LHC and at the Tevatron, it is crucial to obtain 
theoretical predictions for the measured observables which are as accurate as possible. The 
total pair-production cross section has been known at next-to-leading order (NLO) for over 
two decades [15-18]. Some time later also differential distributions [19-21] and the forward- 
backward asymmetry [22,23] were calculated to the same accuracy. As the NLO calculations 
suffer from uncertainties larger than 10%, it would be desirable to extend them to next-to- 
next-to-leading order (NNLO) in a fixed-order expansion in the strong coupling constant. 
Full NNLO predictions require the calculation of two sets of corrections: i) virtual corrections, 
which can be split into genuine two-loop diagrams [24-29] and one-loop interference terms [30- 
32] ; ii ) real radiation, which involves one-loop diagrams with the emission of one extra parton 
in the final state, and tree- level diagrams with two extra partons in the final state [33-36]. 
In spite of progress made by several groups on different aspects of the NNLO calculations 
in the last few years, especially in developing a new subtraction scheme and calculating the 
contributions from double real radiation [37, 38] , a significant amount of work is still required 
to assemble all the elements. 

Another way to improve on the fixed-order NLO calculation (and also the NNLO one, 
upon its completion) is to supplement it with threshold resummation [39,40]. More precisely, 
one identifies a threshold parameter which vanishes in the limit where real gluon emission is 
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soft, expands the result to leading power in this parameter, and uses renormalization-group 
(RG) methods to resum logarithmic corrections in this parameter to all orders in the strong 
coupling constant. For the total cross section, one such approach is to work in the threshold 
limit (3 = a/1 — 4m^ / s — > 0, where s is the partonic center-of-mass energy squared, and f3 
is approximately the velocity of the top (or anti-top) quark. In this production threshold 
limit the top quarks are produced nearly at rest and there are logarithmic terms of the form 
a™ ln m (3 (with m < 2n), which can be resummed to all orders. This has been done to leading- 
logarithmic (LL) order [41-46], next-to- leading- logarithmic (NLL) order [47], and next-to- 
next-to-leading-logarithmic (NNLL) order [48-53]. Note that in this case not only logarithmic 
corrections, but also Coulomb corrections involving inverse powers of (3 occur. An approach 
to all-orders resummation of Coulomb terms can be found in [54] . 

A drawback of this method is that it performs resummation of terms that become important 
in the region j3 — > 0, which however gives a very small contribution to the total cross section. 
At the Tevatron and LHC, typical values of f3 are in the range between 0.4 and 0.9 [55]. 
An alternative approach is to perform the threshold expansion and resummation at the level 
of differential distributions, which are interesting in their own right, and to obtain the total 
cross section by integrating the results. One such method works with the top-pair invariant- 
mass distribution da/dM, where M is the invariant mass of the ti pair. The threshold limit 
for this case of pair-invariant-mass (PIM) kinematics is defined as z = M 2 /s — > 1, and the 
corresponding threshold logarithms are of the form a™ [ln m (l — z) / (1 — z)] + (with m < 2n — l). 
In this limit only soft gluons can be emitted, but (3 is a generic 0(1) parameter, and the top- 
quark velocity need not be small. Systematic resummation and fixed-order expansions of these 
logarithms has been studied in Mellin moment space at NLL order [56-63] and recently also in 
momentum space at NNLL order [55,64], using techniques from soft-collinear effective theory 
(SCET). In general, it can be imagined that the approach based on differential distributions 
captures more contributions than the approach based on the small-/3 expansion, and therefore 
gives more reliable predictions for the total cross section. In [55], we have argued that this 
is indeed the case. We also pointed out that in the effective-theory calculation the argument 
of the soft logarithms involves the ratio 2E s /fi, with 2E S = M(l — z)/y/z the energy of the 
soft radiation. This allowed us to identify a set of power-suppressed terms proportional to 
In z/(l — z), and a numerical analysis showed that incorporating these terms in the fixed- 
order threshold expansion greatly improves the agreement between the approximated results 
in PIM kinematics and the exact results, as was indeed expected based on previous studies of 
Drell-Yan [65] and Higgs production [66,67] in SCET. 

In addition to the invariant-mass distribution of the top-quark pair, the transverse-momen- 
tum and rapidity distributions of the top quark (or anti-top quark) are also interesting. In 
the case of distributions of the top quark, one collects the anti-top quark and extra radiation 
into an inclusive hadronic state X[t\ with total momentum px, and defines the threshold limit 
as p x — > mf. In this limit of single-particle inclusive (1PI) kinematics, only soft radiation 
is allowed, but as in PIM kinematics the parameter (3 is an 0(1) quantity. Soft gluon re- 
summation for this case has been developed in Mellin moment space [68] and applied to the 
top-quark transverse-momentum distribution at NLL order [58] and recently also at NNLL 
order [69], in the form of approximate NNLO predictions. Starting from these distributions, it 
is possible to obtain the total cross section by integrating over the kinematic variables. This 



2 



approach has been taken in [59,60,63,69], and a large discrepancy from the results based on 
PIM kinematics was found. It was argued in [60,69] that the results based on PIM kinematics 
largely underestimate the cross section, while those based on 1PI kinematics more reasonably 
account for the higher-order corrections. In principle, in the threshold limit these two kine- 
matics encode the same soft gluon physics. Any differences between the two cases are due 
to power-suppressed corrections. At realistic collider energies, however, subleading terms are 
in general non-negligible, and one should study them carefully before drawing any definite 
conclusions. In light of the improved behavior of the PIM expansion upon the inclusion of 
the \nz/(l — z) terms observed in [55], it is natural to ask whether including the analogous 
set of terms in 1PI kinematics has a similar effect, and if so, whether it can account for the 
discrepancies between the two types of kinematics observed in previous work. 

In this paper, we extend our effective-field theory approach for PIM kinematics to the case 
of 1PI distributions, performing an NLO+NNLL resummation directly in momentum space, 
and also deriving approximate NNLO formulas equivalent to those from [69]. We then apply 
these results to the transverse-momentum and rapidity distributions of the top quark, as well 
as the forward-backward asymmetry at the Tevatron. As in the case of PIM kinematics, in 
the effective-theory formulation the soft threshold logarithms involve the ratio 2E s //j,. This 
allows us to identify a set of power-suppressed terms proportional to ln(l + s^/m^/s^, where 
s 4 is defined in (4) below and approaches zero in the threshold limit. Through detailed studies 
of differential distributions and the total cross section, we show that although the power- 
suppressed terms in 1PI kinematics can be significant, including the set of power corrections 
identified through the effective field-theory techniques greatly improves the behavior of the 
1PI threshold expansion in regions where the differential cross section is large. In fact, we find 
that the results for the total cross section at NLO+NNLL and approximate NNLO are very 
much compatible within the two types of kinematics, as long as the extra terms appearing 
within the effective-theory analyses are included. 

The organization of the paper is as follows. In Section 2, we review the kinematics of 
top-quark pair production in hadronic collisions and define the threshold regions in 1PI and 
PIM kinematics. Section 3 contains the main new theoretical results of our work. In that 
section we explain how to extend the SCET formalism to 1PI kinematics and assemble all 
of the perturbative ingredients needed for NLO+NNLL and approximate NNLO predictions. 
This includes a calculation of the soft matching function at NLO in a s , which is a necessary 
ingredient for NNLL resummation but was so far unknown. We also demonstrate how the 
extra subleading terms in the threshold expansion appear naturally in our formalism. In Sec- 
tion 4, we define the scheme needed to turn these theoretical results into numerical predictions. 
This includes not only a review of input parameters such at the top-quark mass and parton 
distribution functions (PDFs), but also a definition of the scheme used in evaluating the for- 
mulas at threshold. We then turn to numerical results in Section 5, performing a very detailed 
comparison of the 1PI results with exact results at NLO for the differential distributions, and 
also with the PIM results at NLO and NNLO for the total cross section. Finally, in Section 6 
we give concrete results for the transverse-momentum and rapidity distributions, the forward- 
backward asymmetry at the Tevatron, and the total cross section, both at NLO+NNLL and 
approximate NNLO. We conclude with a summary of our findings in Section 7. 
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2 Kinematics 



We consider the scattering process 

iVi(Pi) + N 2 (P 2 ) t(p 3 ) + t(p 4 ) + X, (1) 

where Aq and N 2 indicate the incoming protons (LHC) or proton and anti-proton (Tevatron), 
while X represents an inclusive hadronic final state. In the Born approximation, two different 
production channels contribute to the scattering process (1): the quark- ant i- quark annihilation 
and gluon-gluon fusion channels. At higher orders in the strong coupling constant, in addition 
to these two channels, there are also contributions from other partonic channels such as quark- 
gluon scattering. These additional channels are power suppressed in the partonic threshold 
region which we will define below, therefore we will not discuss them in this work. Nevertheless 
we will include them in the numerical results when matching with fixed-order calculations. 
The partonic processes which we will analyze in detail are thus 

q(Pi) + qfa) -> t(p 3 ) + t(p 4 ) + X(k) , 

g(pi) + g{P2) t( P3 ) + i( P4 ) + x(k) . (2) 

Note that the hadronic state X in the above equations is different than the state X in (1): X 
contains only the products of the hard-scattering, while X contains also the beam remnants 
from the initial hadrons. The relations between the hadronic and partonic momenta are 
Pi = X\Pi and P2 = x 2 P 2 . At the hadronic level, we define the kinematic invariants as 

s=(P 1 + P 2 ) 2 , h = (P ± - p 3 ) 2 - m 2 t , Ul = (P 2 - p 3 f - m\ , (3) 

while the corresponding quantities at the partonic level are given by 

S = X±X 2 S , t\ = X\t\ , Ui = X 2 U\ , 

M 2 = (p 3 + p 4 ) 2 , s 4 = s + h + «! = (p 4 + kf - m 2 t . (4) 

Momentum conservation implies that at Born level (for k = 0) we have s = M 2 and s 4 = 0. 

Given the kinematic variables above, it is possible to define different threshold regions 
depending on the observables of interest. For example, in the case of the invariant-mass 
distribution of the top-quark pair, it is natural to define the threshold region as s — >■ M 2 . This 
case, which is conventionally called pair invariant- mass (PIM) kinematics, has been analyzed 
in [55,57,59,60,63]. On the other hand, in the case of single-particle inclusive (1PI) observables 
such as the transverse-momentum or rapidity distribution of a single top quark, it is natural 
to define the threshold region as s 4 — > 0. This threshold limit was first studied in [68] and 
is the subject of this work. It should be stressed that, in both the PIM and 1PI threshold 
regions, the top and anti-top quarks are not forced to be nearly at rest. Instead, the top-quark 
velocity f3 = a/I — Am 2 / s is considered a generic 0(1) parameter, neither close to zero nor 
close to unity, which is in fact the region of phase space where the differential cross section 
is largest [55]. This differs from the situation encountered in many calculations of the total 
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top-quark pair production cross section, which are carried out in the limit of vanishing top- 
quark velocity j3 — > [47-49,58-60,63]. We will refer to the kinematic region where (3 — > as 
the production threshold limit. Note that at very large transverse momentum the top-quark 
velocity /3 approaches 1. To deal with such highly boosted top quarks would require a different 
type of effective theory, which has been developed in [70] for e + e~ collisions, but has not yet 
been extended to hadron colliders. 

The common link between these different threshold limits is that they in one way or 
another force the emitted real radiation to be soft. One uses this fact to define a partonic 
expansion parameter which vanishes in the exact threshold limit, where the energy of the 
soft radiation goes to zero, and calculates the partonic hard-scattering kernels as a threshold 
expansion in that parameter. However, unless restrictive kinematic cuts are applied, the 
convolution integral between the partonic hard-scattering kernels and PDFs is sensitive to 
kinematic regions where the partonic threshold parameter is not small. In that case, power 
corrections to the threshold expansion can become important, and these have a different 
structure in each of the three limits mentioned above. In [55] , we analyzed several advantages 
of calculations of the total cross section carried out in the partonic threshold limit for PIM 
kinematics with respect to calculations based on the production threshold limit. Similar 
considerations apply to the partonic threshold limit for 1PI kinematics, so when we discuss 
the total cross section in the present work we focus instead on a comparison between 1PI and 
PIM. 

Several phenomenologically interesting observables are related to the detection of a single 
particle rather than a pair, and 1PI kinematics is suitable to describe such observables. In 
this paper we focus our attention on the top-quark transverse-momentum and rapidity dis- 
tributions, as well as on the total cross section. The double-differential cross section in the 
transverse momentum p T and the rapidity y can be written in the factorized form 



da Wnp-j 



dprdy 3s 



^ dXi ^ ^ f. /Ni ( Xl ^ f } f.^^^^Q^S^^k^^mt^f) , (5) 
j ■ Jxf ln X l Jxf™ X 2 



where the fim are universal non-perturbative PDFs for the parton i in the hadron N, and 
the hard-scattering kernels are related to the partonic cross section and can be calculated 
perturbatively as series in the strong coupling constant. We define their expansion coefficients 
according to 



Cij = a] 



ij 1 



(6) 



The hadronic Mandelstam variables are related to pt and y in the laboratory frame via 

ti = — \fsm± e~ y , u\ = — ^m±e y , (7) 

where m± = \Jp\ + mf . Together with (4), the kinematic variables entering can be 
expressed as functions of pr, y, %i and xi- The lower limits of integrations in (5) are 

-Ui min -Xiti 



mm _ ^ mm _ ' 1 ' ' 



ti L X\S + U\ 
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Finally, the total cross section can be calculated by integrating the double-differential distri- 
bution over the ranges 



1, 1 + Jl - Am 2 s 

< \y < - In V x/ : 

"''"2 l-y/l- 4m 2 Js 



V<PT<^- A 



(8) 



In this paper we will mainly discuss the distributions of the top quark, but our results can 
also be applied to the transverse momentum px and rapidity y distributions of the anti-top 
quark after appropriate replacements. At the Tevatron, charge-conjugation invariance of the 
strong interactions implies that within QCD we have the simple relation 

da da 



(9) 

dp T dy dp T dy PT ^p T ,y=->~y' 

which we will use later on in interpreting the charge asymmetry in terms of a forward-backward 
asymmetry. 

In the threshold limit 54—7-0 only the emission of soft radiation is allowed. In this limit 
calculations are greatly simplified, since one effectively deals with a two-body final state. The 
hard-scattering kernels can be factorized into a hard function and a soft function S^- as 

Cj J (s 4 ,s,fi,it 1 ,m t ,/x / ) = Tr [H" ij -(s / ,*' 1 , u[, m t ,fx f ) Sy(s 4 , s' ,i[,u[,m t , fx f )] + 0(s 4 ) . (10) 

The notation above is meant to emphasize that there are ambiguities in the choice of s', t[ 
and Hi, which can in general differ from the exact Mandelstam variables s, t\ and ui by power 
corrections that vanish at s 4 = 0. For instance, given an explicit result for the hard and soft 
functions, one can always rewrite it using s' + t[ + u[ = or s' + t[ + u[ = s 4 . Although 
the difference is suppressed by positive powers of s 4 , the two choices give different numerical 
results upon integration. In Section 4, we will explain in detail our method for dealing with 
this ambiguity. 

The use of boldface in (10) indicates that the hard and soft functions are matrices in color 
space. The hard function originates from virtual corrections and is the same as that in the 
PIM case [55], once the variables s', t[, u[ are expressed in terms of M and cos#. The soft 
function captures contributions arising from the emission of soft real radiation. It depends on 
the details of the phase-space integrals and is different in 1PI and PIM kinematics. The 1PI 
soft function contains singular distributions in s 4 , which are of the form 



Pu{sa) = 

where the plus-distributions are defined by 



ml 



s 4 



+ 



dsA 







iln»it 



s 4 



m 



g(s 4 



With this definition 
1 



^0 



In" 



s 4 



2 

mi 



+ 



f 



ds 4 -\n n (^)[g(s 4 )-g(0)]. 



s 4 



m 



;n) 



(12) 



ds 4 — ln r 

s 4 



2 

mi 



[g(s 4 ) - g(0)] + 



9(0) 
n+ 1 



In 



n+l 



mi 



(13) 
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Near threshold, these singular distributions lead to a bad convergence of the perturbation 
series. More generally, they give rise to the dominant higher-order corrections to the hadronic 
cross section if the product of PDFs in (5) falls off very quickly away from values of x±,x 2 
where s 4 — > 0. In either case, resumming such terms to all orders can lead to improved 
theoretical predictions. This is the topic of the next section. 



3 Factorization and resummation in SCET 

In [55], factorization and resummation for the case of PIM kinematics was studied in detail. 
Most of that discussion can be carried over directly to the case of 1PI kinematics. In fact, the 
differences between 1PI and PIM kinematics arise solely from the structure of real emission 
in the two cases and therefore affect only the soft function, which must be modified in two 
essential ways. First, the phase-space integrals for real emission in the soft limit change, so 
the fixed-order expansion of the soft function is different from its PIM counterpart and must 
be calculated from scratch. Second, the RG equation for the soft function, derived using the 
RG invariance of the cross section along with the evolution equations of the hard function and 
PDFs, also differs slightly from its expression in PIM kinematics. In what follows we focus 
on how to deal with these two differences with respect to PIM kinematics and otherwise just 
quote results from [55] for the pieces which remain unchanged. In particular, we derive the 
explicit one-loop soft function and the RG equation needed for NNLL resummation in 1PI 
kinematics, and present results for approximate NNLO formulas in fixed-order perturbation 
theory. We also discuss the structure of power corrections to the leading-order term in the 
threshold expansion and explain how a certain set of subleading corrections in S4 appears 
naturally within the SCET formalism. 



3.1 Soft function in 1PI kinematics 

In general, the soft function is related to the vacuum expectation value of a soft Wilson-loop 
operator. To calculate it explicitly, we first generalize the derivation of the differential cross 
section at partonic threshold given for PIM kinematics in [55] to the 1PI case. 

In the limit where extra gluon radiation is soft, the differential cross section can be factor- 
ized as 1 

1 d3 P 3 f d ^P 4 f ^ x e i(pi+P2-P3-P*)-x 



r 1 <m r gft r 

" " 2s (2tt)32£ 3 J {2nf2Ej 



x ^^Tr[H(s\i[,u[ : m u fjL f )W(x : s\i[ : u[,m t: pL f )] , (14) 

where W is the expectation value of the Wilson-loop operator in position space. Since the 
integrand depends on p 4 only through E A = a/I^I 2 + mf and in the exponent, its calculation 
is simplified by going to the rest frame of the inclusive final state t + X in (2), which consists 



lr The functions H and W are summed over the channel indices ij. In order to keep the notation as simple 
as possible, in the remainder of this section we suppress the sum and indices. 
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of the anti-top-quark plus additional soft radiation. In this frame |p 4 | = 0(s 4 /m t ), and we 
can drop the dependence on it in £4 ~ m t . The integral over p 4 then produces a factor of 
(27r) 3 (5^(x), and the exponent depends only on E s = (pi +p 2 — p 3 — p 4 )° = s 4 /(2a/s4 + m 2 ), 
which is the energy of the soft radiation in this particular frame. Using the (^-function to 
perform the integral over d 3 x, and converting the result to a differential cross section in t 1 and 
Mi, we obtain 

da 8n 1 f dx° ( ix°s<i 



i_ r <h^ ( ix°s 4 \ 

n t J 4tt 6XP \2Js 4 + mf) 



dhdui 3s 2 m t J Arc \2^s 4 + m 2 , 

x Ti[H{s',i[,u[,m t ,fi f )W({x°,x = 0),s',i' 1 ,u[,m t , f i f )} . (15) 
We now introduce the momentum-space soft function according to 

Wiuj^s'J'^u'^mt,^) = J ^exp (^r^J W((x°, x = 0), s', i[, u[, m t , fi) . (16) 

The soft function entering the factorization formula is then given by 

S(s 4 ,s',*i,«i,m t ,fi) = — W\—= Sa = ,s'j[,u[,m t ,ii) 



(17) 



Using the results above, one recovers the factorization formula (10) for the hard-scattering 
kernels by noting that 

3s 2 da - 

C i:j (s 4 , s, £1, Mi, m t , «/) = 2 ji ■ ( 18 ) 

07T dtidui 

It is instructive to compare this derivation with that given in [55] for PIM kinematics. 
In both cases, the calculations are simplified by working in a frame where the soft function 
depends only on the time component x°, or in momentum space, on the energy of the extra soft 
radiation (although of course the final results are Lorentz invariant and do not depend on the 
frame). In PIM kinematics this is the partonic center-of-mass frame, while in 1PI kinematics 
it is the center-of-mass frame of the unobserved partonic final state, consisting of the anti-top- 
quark plus additional soft radiation. The difference between these two cases comes from which 
combinations of the momenta are counted as "small" : in PIM kinematics it is Q93 + p 4 ) • k 
against (pz+pi) 2 = M 2 , while in 1PI kinematics it is p 4 - k against p\ = m 2 t . This difference has 
important implications for the structure of power corrections in the two types of kinematics. 
Such power corrections come both from time-ordered products involving the subleading SCET 
Lagrangian and operators, which before phase-space integrations are the same in both cases, 
and from the approximations in the phase-space integrals, e.g. E4 = a/Ip^I 2 + m t ~ m t m 
1PI and E 3 + E 4 = a/|p*3 + p 4 | 2 + M 2 ~ M in PIM kinematics. The corrections from the 
latter source can be quite different: for 1PI kinematics they involve the expansion parameter 
\P4\/mt — s 4 /(2m t A/m 2 + s 4 ), while for PIM kinematics they involve the expansion parameter 
\P3 + PA/M = - z)/(2M) = (1 - z)/(2y/z). 
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A complete analysis of power corrections is beyond the scope of this paper. However, in the 
case of Drell-Yan [65] and Higgs production [66,67] near threshold it was found that by keeping 
the exact dependence u = 2E g in the SCET soft functions analogous to (17) one can reproduce 
a set of logarithmic power corrections involving lnz/(l — z), which are indeed present in the 
analytic results for the fixed-order expansions of the hard-scattering kernels. Keeping such 
terms improved agreement of the threshold-expanded hard-scattering kernels with the exact 
results in QCD. In [55], we followed this procedure in PIM kinematics, using the exact form 
u = 2E g = M(l — z)/y/z instead of u — M(l — z) in the first argument of the soft function. 
Results in this "PIM S cet" scheme include the same type of logarithmic corrections found 
in Drell-Yan and Higgs production, and the numerical results for the threshold expansion 
at NLO in this scheme are indeed improved significantly compared to the traditional PIM 
scheme [59,60,63], which does not include such corrections. In 1PI kinematics, the equivalent 
procedure is to use u = 2E g = s 4 / \/m 2 + s 4 , as we have already indicated explicitly in 
(17). We will refer to numerical results obtained with this choice as being calculated in the 
"IPIscet" scheme. Since the factorization formula (15) is derived in the limit s 4 <C m 2 , it 
would be equally valid to use u = s^jm t in the first argument of the soft function (17). This 
is in fact the choice that has been made in previous calculations in 1PI kinematics [68], and 
later on we will refer to this as the "1PI" scheme. When expanded in fixed-order perturbation 
theory, the two schemes differ through terms involving ln(l + s 4 /m^)/s 4 , which are power 
suppressed in the limit s 4 — >■ 0. However, in our analysis in Sections 5 and 6 we will see that 
these power-suppressed effects can be numerically important, and that the agreement with the 
exact numerical results at NLO is improved in the IPIscet scheme. Furthermore, although in 
this case we do not have explicit analytic results to compare with, we note that the logarithms 
of ln(l + s 4 /raf) appear naturally in the fixed-order NLO calculations of hard gluon corrections 
through terms of the form ln(2Eg//i), see for instance Eqs. (4.16) and (4.17) of [16]. 

We now present the calculation of the soft function in 1PI kinematics at one- loop order. 
The results can be written as 



W^e,**,?) = J2 (19) 



where the matrices are related to products of color generators and can be found in [55]. 
The integrals X[- are defined as 

2^.( Cj x ,|i) = -Mf Vi-vj / d*k 6 (2*)6(k?)8(k°). (20) 

Evaluating these objects in the rest frame of the heavy anti-top quark, we find for the non- 
vanishing integrals 



T - 



2 2/ , s'm 2 \ / , s'm 2 t \ 2 n 2 ( s'm 2 

+ _ Lo _ln * + Lo-ln^ +- + 2Li 2 



^ 3 = - + 2L -^i±^ln* s , 
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T 



+ 2L + 4, 



x 14 — x 24 



T' 
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x 23 



T' 

x 34 



1 1 

^7 + - 



+ 2Li 2 



X( 3 $ U[) , 



7T 

12' 



(21) 



2 In 



2 In 
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7T 
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+ 2Li 2 



2A 



2 7T 2 ' 

— lnx s — 2Lq \xix s + 2 In 2 x s — 4 lnx s ln(l -ij) — 2Li 2 (x 2 ) + — 

6 O 



where (3 t = y/l - 4m 2 t /s>, x s = (1 - &)/(! + A) and 

/i 2 (:r ) 2 e 2 ^ 



In - 



(22) 



The renormalized function is obtained by subtracting the 1/e™ poles in the bare function. 
When performing resummation it is more convenient to introduce the Laplace transform of 
this object, which is defined as 



PCX) 

?(L, s', u\, m t: n) — I duexp 

Jo 



UJ 



= W\\x 



2i 



W(w,s' ',*i,ui,m t ,Ai) , 



e~< E {ie L / 2 ' 



x = ),s ,t 1 ,u 1 ,m t ,fi . 



(23) 



From the second line it is clear that the Laplace-transformed function s is determined directly 
from the position-space soft function through the replacement L — > —L [65]. 

The above expression for the one-loop soft function is new, but we have been able to 
perform two important checks. First, we have verified that our results are consistent with 
the results for real emission in the soft limit given in [16] for the gg channel and in [17] for 
the qq channel. Second, we have made sure that the divergence structure of the one-loop soft 
function is consistent with the RG equation derived in the following section. 



3.2 RG equations and resummation 

The physical cross section should be independent of the factorization scale, which implies 







r 1 dxi r 1 

Jj.min X\ J x n 



d,X2 

x 2 



din/// 



[fi/Ni{x U Hf) fj/N a (x 2 ,flf)] 



Cij(s 4 , s,ti,iti,m t ,///) 
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+ fi/Ni (Xl, /If) fj/N 2 (x2, /if) ^~^ C '^( g4 ' S ' *1' "1' m *' Vf) } • ( 24 ) 

Compared to the PIM case, the terms arising from the derivatives acting on the PDFs are 
slightly different. To understand their structure, first consider the term where the derivative 
acts on fj/Ni(x2,fJ>f)- Under the dynamical assumption of steeply falling PDFs, the DGLAP 
evolution equations can be simplified by keeping only the leading terms for x — > 1 in the 
Altarelli-Parisi splitting functions, in which case the evolution becomes diagonal and reads 



d f 1 dx 

fj/N 2 (y, n) = j — Pjj (x) f j/N2 (y/ x, n) , (25) 

Jy X 



dlnfi 



with Pjj(x) given by 



2P (a ) 

PjiW = ( r P \ + 27?(«.) 5(1 - x) • (26) 

(1 - X) + 



We then obtain (up to power-suppressed terms) 

• (27) 



1 dxi 

Ji/NAxi, /if) 

min db\ 



— I T" hm( x *lii Vf) P n(0 Cij(s 4 , s, ti, mi, m t , Hf) 

fin X 2 J X2 C, 



To derive the evolution equation for the soft function, we arrange the integrations such that 
the Altarelli-Parisi kernel acts on the CV, rather than on the PDF. After some manipulations, 
the term in the square brackets above can be written as 

\ — fj/N 2 (x 2 ,Vf) / —Pjj^Ci^s'^is^i^iu^rnt^f) (28) 

J^mm X2 Jx 2 nln /X 2 S 

with s' A = ti + £(s + iii). Changing the integration variable from £ to s' A , taking the threshold 
limit s 4 y 0, and using the identity 

(-*l)l(S4-sy/(-fi)]+ Jo S 4 -S 4 

we can convert this term into a form from which the RG equation for the soft function is more 
easily derived. The analogous procedure is then used for the term where the derivative acts 
on the other PDF. Finally, we use the RG equation for the hard function 

-H(s' \i!i,u'i,m u ii) = T H (s\i[,u[, fi) H(s\i\,u\,m t , fi) 



dlnjj, 

+ H(s, i[, u[, m t , n) T ] H (s, i[, u[, m t , //) , (30) 



with 



T H (s, i[,u[, m t , /-0 = r cusp (a s ) ( ln - iir ) 1 + j h (s, i[,u[, m t , a s ) , (31) 
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where r cusp is equal to C_f7cus P for the qq and C^Tcusp fo r the gg channel, and the matrices 
7^ can be found in [55]. 

Assembling the different pieces and enforcing (24), the evolution equation for the momentum- 
space soft function reads 



d 



dlnfi 



S(s 4 ,s w^m^/i) = - 



2r cusp (a s ) In + 7 st (s', i[, u[, m t , a s 



S(s 4 , s ,t 1 ,ii 1 ,m t , n) 



- 5(s 4 ,s / ,t' 1 ,M / 1 ,m t ,/i) 

(«.) / 
Jo 



s 4 



2r cusp (a s ) In — ^- + 7 s (s', t[, u[, m t , a s ) 



-AT 



cuspV 



where we have defined 



, S(s , A ,s , ,i' 1 ,u' 1 ,m t ,(i) - S(s 4 ,s',ti, u'^rrit,^) 
ds A 

S4 - S' A 



"I 2 

s m 



(32) 



Y(s', t' 1: u' 1: mt, a 8 ) = j h (s' : t[, u' 1: m t: a s ) + 2>f{a a ) + T cusp (a s ) In ^ 1 . (33) 

This evolution equation is of the same form as for the PIM case, but the soft anomalous 
dimension is modified by the logarithmic term in (33), which can be traced back to the 
different form of the collinear evolution terms after arranging the integrations as appropriate 
for 1PI kinematics (this extra term vanishes in the production threshold limit ft — >■ 0, where 
PIM and 1PI kinematics agree). Therefore, we can use the expression for the resummed soft 
function derived in [55], taking into account the changes in the anomalous dimension and soft 
matching function, and the fact that now oj = s 4 / \J s 4 + m\ sets the mass scale in the soft 
logarithms. The resummed soft function is then given by 

*7 2 

s m 



-4S(fj, 8 ,fj,f) + 2a r (/i s ,/i/)ln^-^ + 4a 7 «(/z s ,///) 

*iui 



S(s4,s',t' 1 ,u' 1 ,m t ,iJ 1 f) =exp 

x u f (s, i[, u[, mt,Hf, /i s ) s(d v , s', i'^u^ m t , /x s ) u(s, i'^u'^ m u jif, fj, s ) 



x 



s 4 



S4 \^m 2 t + s^ s J T(2r/) 



2?) e -2 lEV 



(34) 



r)=2a r (ti s ,fi f ) 

and combining this with the solution for the hard function, the final result for the resummed 
hard-scattering kernels is 



C(s 4 , s , t ls u ls m t ,n f ) = exp 



x Tr 



U(s', t[, u[,m t , n h , fx s ) H(s, i[, u[,m t , fi h ) U f (s', i'^u'^ m t , n h , /x s ) 



x s(d v ,s',i / 1 ,u' 1 ,m t ,n s ) 



s 4 



s 4 Ky/mt + SiHsJ r(2?7) 



2r? e -2 7E r, 



(35) 
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The result (35) has already been given in [55], where the explicit definitions of the RG expo- 
nents ar and and the evolution factors U and u can be found. The notation is such that 
one must first take the derivatives with respect to rj appearing in the first argument of the 
Laplace-transformed soft function s, defined in (23), and then set r\ = 2ar(// s , ///). For values 
of the scale where 77 < 0, one must use analytic continuation to interpret the formula in terms 
of plus-distributions. 

Our result for the resummed hard-scattering kernels is equivalent to the Mellin-space re- 
summation formula from [68] when expanded to any fixed order in a s , if we approximate 
a/ m 2 + s 4 ~ m t . However, the scale choices used in Mellin-space resummation are typically 
such that the one encounters a Landau-pole ambiguity in the evaluation of the all-orders for- 
mula, upon inverting the Mellin transform and integrating over s 4 . One way of dealing with 
this, as done in [59,60,63,69], is to instead use the resummation formula only to construct 
approximate fixed-order expansions at NNLO. Another way is to instead view ji s as a function 
of the observables and y, and choose it in such a way that the perturbative expansion of 
the soft function at fi s is well behaved. With such a choice of /x s , one can still evaluate the 
all-orders resummation formula, but without encountering Landau-pole ambiguities. We will 
describe how to construct the fixed-order expansion to NNLO in the next section, and then 
compare the two methods in more detail in Section 3.4. 



3.3 Approximate NNLO results 



In addition to the resummed formulas derived above, it is sometimes useful to construct 
approximate NNLO formulas in a fixed-order expansion. In our formalism, this is achieved by 
evaluating the expressions 



C(s 4 , s', i[, u[, m t , fj) = c(d v , s, i'^u^, m t , ji) 



S 4 



r (2f7) S4 \^m 2 t +s 4 fi & 



2,1 



jj=0 



where 



c(L, s', i[, u[, m t , jj) = Tr H(s', i[, u[,m t , jj) s(L, s 1 , i'^u'y, m t , jj) 



(36) 



(37) 



Using the ingredients of our NNLL resummation formula, it is possible to use RG methods to 
determine at NNLO the coefficients of all powers of L in c, as well as the /i-dependent part 
of the constant piece. To do so, we use the same methods as in [64]. We then convert the 
derivatives with respect to the auxiliary parameter r\ into distributions in s 4 defined in (11), 
which can be easily done by using the following replacement rules: 



1 

L 

L 2 
L 3 



5(s 4 ) , 

2-Po(s4) - S(s 4 ) L m , 
8Pi(s 4 ) - 4L m P (s 4 ) + S(s 4 ) ( L 2 m - 



2tt 5 



4L 4 



3 ) s 4 

24P 2 (s 4 ) - 2AL m P l {s 4 ) + (QL 2 m - 4vr 2 ) P (s 4 ) + 5(s 4 ) {-L 3 m + 2n 2 L m + I6C3) 
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6L 4 

s 4 



-L 4 + 2 In 



^2 1 



m 2 fi 2 



-> 64P 3 (s 4 ) - 96L m P 2 (s 4 ) + (48L 2 t - 32tt 2 ) P 1 (s 4 ) + (-8L^ + 167r 2 L m + 128C 3 ) P (s 4 ) 

{ 4tt 4 " 
+ 8(s A ) Lt - An* Li - 64( 3 L m + 



15 



8L 4 
s 4 



s 2 
s 4 



<? 2 

2 6 4 



L\ - 3L 4 In + 3 ln z — ^ - 2^ 



m 2 /i 2 



(38) 



where L m = ln(/i 2 /m 2 ) and L 4 = ln(l + s 4 /m 2 ). 

The final result for the hard-scattering kernels at NNLO can be written as 

C (2) (s 4 , s', i[, u\, mt,n) = D 3 P 3 (s 4 ) + D 2 P 2 (s 4 ) + D 1 P^) + D P (s 4 ) + C 5(s 4 ) + R(s 4 ) , 

(39) 

where the coefficients D , . . . , D 3 and C are functions of the variables §', t[, u[, m t and fx. 
The explicit results are quite lengthy and are contained in a computer program which can be 
downloaded together with the arXiv version of this paper. The analytic expressions for Di 
(i = 1, 2, 3) were first derived in [60] starting from resummed formulas in Mellin moment space. 
We have compared with those results and found agreement. In a recent paper the coefficient 
D was also determined [69] , but its explicit form was not reported there. We are thus unable 
to compare this term. The regular piece R(s 4 ) collects terms involving L 4 , which arise from 
choosing uj = s 4 / \/ m 2 + s 4 in the argument of the soft function. As noted in Section 3.1, 
dropping R{s 4 ) recovers the 1PI scheme used in earlier work, for instance [69]. The Co term, 
on the other hand, is ambiguous since only its scale-dependent part is exactly determined, and 
one needs to specify which contributions are included there. One contribution to C comes 
from the conversion of powers of L in according to (38), which is determined exactly, 
therefore the ambiguity comes from the constant term of c^ 2 \ which is 



g(2), 



(0) = Tr \h<V s«(0) + s(2)( ) + s^(0) 



(40) 



where we have suppressed the dependence on §', i[, u[, m t and /i for convenience. In the three 
terms above, the first term is known exactly and is therefore included, while the constant term 
of the two-loop soft function in the second term is unknown. As for the two-loop hard function 
in the third term, one can determine its scale-dependent part and include it in the formula, as 
was done for the PIM case in [55]. In that case there is a freedom to write the result in terms 
of ln(/i //i), where the scale /i can be chosen as or m t . However, by including these extra 
/^-dependent terms one runs the risk of artificially reducing the scale dependence, rendering 
it an ineffective means of estimating theoretical uncertainties. Therefore, we will not include 
these terms here and instead drop the contributions of the two-loop hard function. Later, 
when we compare the numerical results for the total cross section derived using PIM and 1PI 
kinematics, the PIM numbers will also be computed in the equivalent way, i.e., by dropping 
the two-loop hard function completely. 
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3.4 Resummation vs. NNLO expansions 



In our numerical studies later on we will typically give results from both resummed pertur- 
bation theory and the approximate NNLO formulas. Although the perturbative information 
used in these formulas is the same - the NLO matching functions and the NNLO anomalous 
dimensions - the implementation and philosophy is different. 

In resummed perturbation theory, one views the soft function as depending on the distri- 
butions 



as well as £(54). The logarithmic corrections are in general considered large compared to the 
5-function term (and of course subleading terms in s 4 ), but it is assumed that with a proper 
choice of fj, s they can be treated on the same footing, so that the soft function at this scale can 
be reliably calculate d in fixe d-order perturbation theory. This would obviously be the case if 
we chose fi s = s 4 / a/s 4 + mf , since then the logarithmic corrections would vanish, but in that 
case the running coupling in the soft corrections would, for some value of s 4 , be evaluated 
at the Landau pole, which would spoil the clear separation between perturbative and non- 
perturbative physics accomplished by using the effective field-theory formalism. Therefore, as 
usually done in momentum-space resummation [55,65-67,71,72], we will view the soft scale as 
a function of the observables pt and y, and choose it based on the convergence of the physical 
cross section. In particular, we study corrections from the soft function to the cross section as 
a function of fi s , and choose the numerical value of the scale as the point where the correction 
is minimized. The formulas then sum logarithms of the numerical ratio /z s //i/, where fi s is the 
dynamically generated soft scale. The same reasoning applies to the choice of the hard and 
factorization scales, and an advantage of the resummation formalism is that the three scales 
can be varied independently as a way of estimating perturbative uncertainties. 

This approach should be contrasted with that based on approximate NNLO formulas. 
In using such an approximation, one assumes that the logarithmic corrections from the P' n 
distributions account for the bulk of the NNLO corrections in fixed order at an arbitrary 
factorization scale /i/, and at the same time that the corrections at NNNLO and beyond, 
evaluated at that scale, are small enough that the perturbative series is well behaved. From this 
point of view, the resummation formalism is just a useful tool for constructing the approximate 
fixed-order expansion, and no physical significance is given to the soft scale fj, s , on which the 
final answer does not depend. 

It is worth emphasizing that the fixed-order expansion of the NLO+NNLL formulas to 
NNLO in a s is not exactly equivalent to the approximate NNLO formulas from the previous 
section. The direct expansion of the NLO+NNLL formulas to NNLO contains explicit de- 
pendence on the scale /i s and also a different pattern of plus-distributions compared to the 
approximate NNLO formulas. For instance, the approximate NNLO formula contains P 3 dis- 
tributions, but the expansion of the NLO+NNLL formula to NNLO in fixed order contains 
at most P2 distributions, as required by the Altarelli-Parisi equations. This aspect of the 
calculation is discussed in greater detail in the Appendix. 




(41) 



15 



The optimal method for including the higher-order perturbative effects is not entirely 
clear without further information. The reason is that it is not possible to tell whether the 
logarithmic corrections from the partonic threshold region can be considered large, or how 
to minimize them with a proper choice of /i s , until after the integration over s 4 . After that 
integration, they give large perturbative corrections to the differential cross section if the PDFs 
fall off very quickly away from the region where s 4 — > 0, an effect referred to as "dynamical 
threshold enhancement" [65]. Since the PDFs are not known analytically, it is only possible 
to assess the extent to which the corrections in the partonic threshold region are dynamically 
enhanced through a numerical study. We thus return to a discussion of the relative merits 
of resummation vs. approximate NNLO expansions after our numerical studies in Sections 5 
and 6. 

4 Numerical implementation 

We now describe the numerical implementation of our results. A central issue is that there are 
power-suppressed ambiguities in the evaluation of (5) and (10) away from the exact threshold 
limit, and we describe in detail how we deal with these ambiguities in the following. We also 
use this section to summarize some of the inputs used throughout the rest of the paper, in 
particular the PDF sets, the strong coupling constant, and the top-quark mass. 

As pointed out in Section 2, there are power-suppressed ambiguities in the choice of the 
variables s', t[ and u[ of the hard and soft functions. Apart from when it appears in the 
(^-function or plus-distributions, in the perturbative calculation of the hard and soft functions 
one can set S4 = everywhere and use s' + 1[ + u\ — to rewrite the hard-scattering kernels 
in many different forms. While these are all formally equivalent in the threshold limit S4 — > 0, 
they change the functional dependence of the hard-scattering kernels on x\ and x 2 , so the 
integration in (5) gives different results for the pieces multiplying plus-distributions in s 4 . 
Moreover, since one typically trades either x\ or x 2 in favor of S4 as an integration variable, 
another obvious choice is to use s'+t^+u^ = S4 before integration, again leading to numerically 
different answers which are nonetheless equivalent in the threshold limit s 4 — > 0. 

Our method of fixing this ambiguity is as follows. First, we enforce s' + 1[ + u[ = in the 
hard-scattering kernels, and use this to eliminate either t[ or u[ as an independent variable. 
We then define the two cross sections 



da 1 



167rp T 




■xi(s+*i)+tti 



dp T dy 



S 4 — X\ti 



x fi/ Nl (xufif) fj/N 2 (x 2 (s4),fx f ) C i:j (s 4 , s', i[, -s - i^mt, ///) , 



(42) 



da 



16np T 




•X2{s+U 1 )+t 1 



ds 4 



dp T dy 



s 4 - X 2 Ui 



x fi/NAM-^),^/) fj/N 2 {x 2 ,Hf) Cij(s A) s', -s' - u' 1: u' 1: m t: iJ, f ) . 



(43) 
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We have changed variables from x 2 or Xi to s 4 in the two equations, respectively, so that 

S 4 — X2U1 , , S4 — X\t\ 



Xi(s 4 ) 



x 2 (s 4 ) 



%2-S + t\ ' " X\S + Mi 

Finally, we drop all dependence on s 4 in the hard-scattering kernels by using 

t[ = i\ = xiti , s' = XiX2(0)s 

in (42), and 

l[ — ill — X2U1 , s = xi(0)x2S 



(44) 



(45) 



(46) 



in (43). It is easy to see that with this choice a 1 and a u are not necessarily the same, although 
the difference is power suppressed. We shall take the average of the two as the final result for 
the differential cross section: 



da 



dp T dy 



do 1 da u 
+ 



dp T dy dp T dy 



(47) 



In this way, the rapidity distribution in the gluon channel is invariant under y — > —y, as it 
should be, and the relation (9) is preserved. 

The scheme above specifies our procedure for the numerical evaluation of the threshold 
formulas for the differential distribution in px and y. In the next section we will also study 
the differential distribution in (3. Rather than define yet another scheme, we will calculate 
this through an exact change of variables and integration orders in (47). The differential cross 
section in /3 then takes the form 



da 4ir/3 



d(3 3sm 



1 

X 2 



+ 



E 



1,] 



/-s(l-/3)/2 ^ ns+k+sml/h 
dii / ds 4 Ci j (s4,s' t ,i' 1 ,-s' t -i' 1 ,m t ,fi f ) 
■s(l+/3)/2 JO 



-S(l-/3)/2 



du\ 



-s(l+P)/2 

where s = 4m|/(l — /3 2 ), and 



s+iii+sm^/ui 



ds 4 Cij(s 4 , s' u , -s' u - u[, u[, m t , fi f ) 



(48) 



si — s 



in the first term in the bracket, while 



s 4 - ti 



-u 4 

S 4 — U\ 



i[ = h 



(49) 



u 1 = Ui 



(50) 
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in the second. We have also defined the parton luminosity function 



fij(y,Vf) 



1 dx 



x 



fi/Ni (x, fi f ) f j/N2 (y/ x, n f ) 



(51) 



Since the change of variables is carried out exactly, (47) and (48) give the same result for the 
total cross section when integrated over. 

We should emphasize again that the procedure above is by no means unique. In [59], for 
instance, the scheme was instead specified at the level of the ft distribution. In particular, 
the hard-scattering kernels as a function of t'^u'^s' were written in a form specified in the 
Appendix of that paper, u[ was eliminated as an integration variable using s' + t[ + u[ = s 4 , 
and then the integration over s 4 was carried out followed by that over i^. 2 We will discuss the 
numerical differences which result from using such an equivalent procedure when we compare 
with previous literature in Section 6.4. 

Having specified our procedure for evaluating the formulas in the threshold region, we 
next clarify how to match the results with fixed-order perturbation theory at NLO. The exact 
results contain the perturbative corrections to our formula which vanish in the limit s 4 — > 0, 
and to obtain solid phenomenological results it is important to include them. In resummed 
perturbation theory, we achieve NLO+NNLL accuracy by evaluating differential cross sections 
according to 



da 



NLO+NNLL 



da 



NNLL 



_|_ ^(jNLO, subleading 



da 



NNLL 



+ da 



NLO 



^ NLO, leading 



where da NLO is the exact result in fixed order, and dcr NLO ' leading | M/ = da 



NNLL I 



(52) 



captures 



the leading singular terms in the threshold limit s 4 — > at NLO. To obtain approximate NNLO 
results in fixed order, we simply compute 



rfa NNLO, approx = ^NLO + ^(2), approx ^ 



(53) 



where da^' a PP rox is the approximate NNLO correction to the differential cross section obtained 
using the coefficient function (39). 

The size of the power corrections contained in parentheses in the second term of (52) 
depends on the region of phase space in which the differential cross sections are evaluated, 
and also on the interplay of the hard-scattering kernels with the PDFs. The power corrections 
are expected to be small when s — > 4m^_, since in that case s 4 — > 0. However, experiments 
do not typically reconstruct s as an observable. For more interesting differential distributions 
the limit s 4 — > can be enforced via a restriction to the machine threshold, for instance by 
requiring that m_|_ — > y/s/2 for the pt spectrum, but in this case the differential cross section 
would be extremely small. Away from such special kinematic regions the leading terms in 



2 We are grateful to Sven Moch for the clarification of this point, and for providing the numerical code used 
in the comparison in Section 6.4. 
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Order 


PDF set 


ot s [M z ) 


LO 

NLO, NLL 

NNLO approx, NLO+NNLL 


MSTW2008LO 
MSTW2008NLO 
MSTW2008NNLO 


0.139 
0.120 
0.117 



Table 1: Order of the PDFs [74] and the corresponding values of the strong coupling used for 
the different perturbative approximations. 

the partonic threshold limit are dominant only if a dynamical enhancement occurs, because 
the product of PDFs appearing in the cross section falls off sharply away from the region 
where s 4 — > 0. Some systematic studies of the conditions under which threshold enhancement 
is effective have been made in [65,73], but in the end this is largely a numerical question, 
which must be examined on a case-by-case basis. Given its importance, we address the issue 
of threshold enhancement and power corrections to the threshold expansion in some detail in 
the next section. 

We end this section by summarizing some of the inputs needed for numerical predictions, 
which we shall use in the rest of the paper unless otherwise indicated. The results depend 
on the input parameters m t , a s , and the PDFs. We shall use the MSTW2008 PDF sets [74], 
and change the PDF set according to the order at which the perturbative hard-scattering 
kernel is evaluated, as indicated in Table 1. The running couplings a s (/i) are taken in the MS 
scheme with five active flavors, using one-loop running at LO, two-loop running at NLO, and 
three-loop running at NNLO. By default we take m t = 173.1 GeV in the pole scheme. Using 
a more physical quark mass such as that in the MS scheme can lead to a better convergence 
of the perturbative series [50,55], but for the purposes of this work we shall leave that issue 
aside. 

Finally, in the following we do not consider the theoretical uncertainty on the top-pair cross 
section induced by the error on a s (M z ), which is an input parameter. However, a recent study 
employing CTEQ PDFs shows that this additional theoretical uncertainty can be as big as 4% 
for the total top-quark pair production cross section at the LHC [75]. Such an uncertainty is 
therefore not negligible in comparison to the residual scale uncertainty in NNLO and NNLL 
calculations, and it will need to be considered when comparing data and theoretical predictions. 
An analysis of the effects of the a s uncertainty on several cross sections of interest at the 
Tevatron and at the LHC was also carried out by the MSTW collaboration [74] . However, the 
top-quark pair production is not among the processes considered in that work. 

5 Threshold enhancement and power corrections 

We now proceed to study threshold enhancement and the numerical importance of power 
corrections to the factorization formula (10). The goal is to examine under which conditions 
the higher-order corrections dominating in the limit S4 — > can be expected to give a good 
approximation to the full result. Given its importance, we approach this question from several 
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different angles. In Section 5.1 we compare the leading terms in the 1PI and IPIscet threshold 
expansions for the pr and rapidity distributions with the exact results at NLO in QCD. We 
then gain more insight into the patterns observed there by studying the differential cross section 
with respect to /3 — y/l — 4mf/s, which is obtained from the parton luminosity functions 
multiplied by the total partonic cross section in Section 5.2. We make general comments on 
the pattern of power corrections in 1PI and PIM kinematics, and also investigate the extent 
by which corrections stemming from logarithmic plus-distributions are enhanced compared to 
other terms. 

5.1 Transverse- momentum and rapidity distributions at NLO 

Our first test of the threshold expansion in 1PI kinematics is to compare its predictions for 
the transverse-momentum and rapidity distributions with the exact ones at NLO in QCD. In 
other words, we study whether the leading singular pieces in the 1PI or 1PI S cet threshold 
expansions provide a good approximation to the exact NLO hard-scattering kernels, so that 
the power corrections contained in the parentheses of the second term in (52) are small. The 
results of this comparison can be found in Figure 1, where the transverse- momentum and 
rapidity distributions at the Tevatron and the LHC with = 7 and 14 TeV are displayed, 
for the choice [if — 2m t . To compute the exact QCD corrections we rely on the Monte Carlo 
programs MCFM [76] and an internal NLO version of MadGraph/MadEvent 3 [77,78], whereas 
to implement the leading pieces of the threshold expansion in IPIscet and 1PI we have used 
the procedure described in Section 4. All inputs, including PDFs and m t , can also be found 
in that section. It is obvious from the figures that the IPIscet approximation does better 
than the 1PI approximation in reproducing the exact QCD results. Since the three curves 
differ only through the NLO corrections to the hard-scattering kernels, we can compare them 
much better if we isolate these pieces. We have done so for the distribution in Figure 2, 
in this case for two different values of [if. We then see that at the Tevatron the IPIscet 
approximation works remarkably well over the full range of [if. At the LHC, the IPIscet 
approximation reasonably reproduces the correction at [if — m t , but does relatively poorly 
in reproducing the correct [if dependence. The 1PI approximation significantly overestimates 
the true result at all three values of [if, both at the Tevatron and at the LHC. 

Given that the numerical differences between the IPIscet, 1PI, and exact results are due 
solely to sub leading terms as s 4 — » 0, we conclude that power corrections to the pure threshold 
expansion can be sizeable at NLO. At the Tevatron, where the qq channel gives the largest 
contributions, the extra terms related to L4 included in IPIscet account for the dominant 
power corrections. They are also important at the LHC, where the gg channel dominates, but 
so are other corrections, which cannot be obtained in our formalism. 

We have focused on the region of pr where the differential cross section is largest. In 
principle, our results can also be used to predict the high-p^ tail of the distribution. In the 
limit m± — > y/s/2, for instance, S4 — > 0, so threshold expansion is bound to work well when 
compared to the exact NLO result. However, in such kinematic regions the differential cross 
section is so small that it is essentially unobservable, and the top-quark is so highly boosted 

3 We are grateful to Rikkert Frederix for providing the code. 
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Figure 1: The transverse-momentum and rapidity distributions at NLO at the Tevatron with 
y/s = 1.96 TeV and at the LHC with = 7 and 14 TeV. 
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Figure 2: The NLO corrections to da/dpx- 
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that m t <C V§, so the power counting in the effective theory would need to be modified. A 
more interesting region would be up to around 400 GeV at the Tevatron, and up to around a 
TeV at the LHC. We will include the higher-p^ region for the Tevatron in the phenomenological 
studies. For the LHC, however, s 4 can be on average rather large at such values of p T , so 
power corrections to the gg and qq channels can become significant, and the qg channel can 
also give non-negligible contributions. Given these problems, we will not study the high-p^ 
distributions at the LHC. 

This simple NLO study is instructive, but in the end the real issue is how well the threshold 
approximation is expected to work at NNLO. On the one hand, at NNLO one encounters 
plus-distributions enhanced by up to three powers of logarithms, so it is not unreasonable to 
expect that the power corrections are of less relative importance than at NLO. On the other 
hand, at NLO the coefficients multiplying both the P n distributions and 5-function terms are 
known exactly, while at NNLO only those multiplying the P n distributions are available. It is 
therefore difficult to anticipate the behavior of the threshold expansion at NNLO based only 
on its behavior at NLO. However, we can gain additional insights through the studies of the 
j3 distribution in the next subsection. 



5.2 The (3 distribution and total cross section 



The differential cross section with respect to /3 — \/l — 4mf / s can be expressed as 



da 
~df3 



8/3 



s(i-/3 2 y 



Am, 



(54) 



The coefficient functions are proportional to the total partonic cross sections, and are 
obtained from our results by comparison with (48). We define expansion coefficients for these 
functions as 



/ y = if + AnaJ^ + (47ra s ) 
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This differential cross section is not measured experimentally, but from the theoretical per- 
spective it is convenient because answers for the partonic cross section are known analytically 



to NLO [18], and for the functions f^' 1 ^ and f-j' z> multiplying the scale-dependent logarithms 
to NNLO [50]. Moreover, this cross section is calculable in both 1PI and PIM kinematics, 
so it gives us a way of directly comparing the threshold expansion and power corrections for 
these two cases. The results must agree in the limit /3 — > 0, since in that case gluon emission 
is soft, but beyond that they receive a different set of power corrections, so the agreement of 
the two approximations with each other is one way of testing whether these power corrections 
are under control. 

We begin with a study of the NLO corrections, similar to that performed for the p T 
distribution in the previous section. In this case, however, we can also include the results from 
PIM kinematics, which we calculate using the expressions in [55]. In Figure 3 we compare 
results for the a s correction to (54) obtained within the different expansions. In contrast 
to the p T spectrum, we compare the quark and gluon channels separately and look at the 
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corrections in a kinematic range from the production threshold at /3 — to the machine 
threshold at /3 max = a/1 — 4m^/s. As anticipated from the results of the previous section, 
the IPIscet approximation works quite well in the qq channel, but somewhat worse in the 
gg channel, especially at the LHC. The 1PI results overestimate the exact corrections in all 
cases. As for the results in PIM kinematics, they are lower than the exact results in all cases, 
but the power-suppressed terms included in PIMscet bring the leading terms in the threshold 
expansion closer to the full result. The PIMscet approximation is slightly worse than IPIscet 
in the qq channel, and slightly better in the gg channel, but the differences are not major. 

As mentioned above, all of the approximations have the same leading-order expansion in the 
limit (3 — > 0. The power-suppressed effects accounting for the differences between the curves 
start to become noticeable at (3 ~ 0.2. After that point, there is more phase space for hard 
gluon emission and the size of the power corrections increases. Evidently, the subleading terms 
included in the IPIscet and PIMscet approximations account for these power-suppressed 
terms in part, although not completely. This is especially noticeable in 1PI kinematics, where 
the power corrections are generically more important than in PIM kinematics, a point we will 
return to below. The power corrections become progressively more important at higher values 
of the collider energy, since then the luminosities are larger at high /3 and the differences in 
the partonic cross sections in that region are magnified. This is most easily seen by comparing 
the LHC results at the two different collider energies, where one observes larger gaps between 
the approximations at 14 TeV than at 7 TeV. A careful examination of the results at the LHC 
with ^/i = 14 GeV also shows a feature not obvious in the other cases: at very high values 
of (3 close to the endpoint, the exact correction in the gg channel remains a positive number. 
In fact, the exact NLO correction to the partonic cross section in the gg channel tends to a 
positive constant at very high s [15], while that for the qq channel, and also the threshold 
approximations to the corrections in PIM and 1PI kinematics in both channels, approach zero 
at high s. This feature is only visible at the highest collider energy, because otherwise the 
high-/? cross section is completely damped by the luminosities. For the same reason, the qg 
channel, for which the partonic cross section also tends to a positive constant at high s, can 
become important at high (3. We show the a s correction from this channel at the LHC in 
Figure 4, for three different values of [if (at the Tevatron, the contribution from this channel 
is still very small.) For higher values of (3, and especially at lower values of /iy, it can be as 
important as the gg channel, even though it is suppressed in the limit s 4 — > 0. We discuss this 
effect at the level of the total cross section in more detail in Section 6. 

In contrast to the transverse-momentum and rapidity distributions, for which little is 
known beyond NLO, for the (3 distribution we can also perform some comparisons with other 
results at NNLO. First of all, results for the (3 distribution at the same level of NNLO approx- 
imation, but in PIM kinematics, were obtained in [55], and the agreement of these results with 
each other gives some information about the size of power-suppressed terms. In addition, the 
NNLO corrections proportional to the scale-dependent logarithmic terms in the f3 spectrum 
(54) are known exactly [50]. This gives us an opportunity to also compare both types of 
kinematics with an exact result beyond NLO. 

The NNLO corrections proportional to the scale-dependent logarithmic terms in the f3 dis- 
tribution (54) within different approximations are shown in Figure 5. The results are obtained 
by dropping the scale-independent coefficient /( 2 '°) from (55). Since these contributions would 
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Figure 3: The a s corrections to da/df3 for the different approximations mentioned in the text, 
with fif = m t . 
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Figure 4: The a s corrections to dcr/df3 from the channel, for /if = m t (solid), /if = m t /2 
(dashed) and /if — 2m t (dotted). 

vanish for /i/ = m t , we have chosen in this case fif = 2m t instead. We have evaluated the 
exact results in QCD using the formulas from [50], and the PIM S cet and PIM results using 
those from [55]. To obtain the threshold expansions for this comparison, we have differed 
from our normal scheme for including the NNLO corrections described at the end of Sec- 
tion 3 by including the /z-dependent pieces of the two-loop hard function, and by rewriting all 
/i-dependent logarithms in the form ln(m^/yU 2 ). Concerning the agreement between the 1PI 
approximations and the exact results, one sees the same qualitative behavior as at NLO. The 
1PI SCET results are consistently a better approximation than the 1PI results, especially at the 
LHC, where the power corrections are large at higher values of (5. As for the PIM results, the 
PIMscet approximation fares slightly better than the PIM results in the gluon channel, but 
slightly worse in the quark channel. In any case, the differences between the PIMscet and 
PIM results are much smaller than those between the IPIscet and 1PI results, which can be 
taken as an indication that the power corrections are smaller in PIM than in 1PI kinematics. 
However, once the extra corrections unique to the 1PI SCET scheme are taken into account, 
the results in these two types of kinematics are very much compatible with one another and 
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provide a good approximation to the exact results. 

The NNLO corrections from the scale-independent pieces /( 2 '°) to the (3 spectrum (54) 
within the different PIM and 1PI approximations are shown in Figure 6, for the choice /if — m t . 
For these pieces it is not possible to make a comparison with an exact result, but we can make 
a couple of comments based on the agreement of the different approximations with each other. 
As before, the difference between the two PIM schemes is small compared to that between 
the two 1PI schemes, indicating that the power corrections in PIM kinematics are smaller, 
and the difference between the 1PI S cet and PIM S cet results are much reduced compared to 
the difference between the 1PI and PIM results. In general, the NNLO corrections in the 1PI 
approximation are much larger than any of the others. 

The explicit results from these studies all point to the fact that 1PI kinematics is more 
susceptible to power-suppressed effects than PIM kinematics. We can gain more insight into 
this observation through a very simple analysis. As discussed in Section 3, the leading power 
corrections in 1PI kinematics are related to the partonic expansion parameter A = E s /m t , 
where 2E S ~ s 4 / a/ + s 4 is the energy of extra soft radiation in the partonic scattering 
process. For the case of PIM kinematics, the equivalent parameter is A = Ef 1M /M, where 
2Ef lM = M(l — z)jyfz. We can quantify in part the relative size of these parameters as a 
function of j3 by evaluating the mean value 




Wipi = / dti ds 4 . / / dt i ds * ( 56 ) 



in 1PI kinematics, where the appropriate integration range can be read off from (48), and the 
analogous expression 



(A) PIM = / d * _ / / dz (57) 




in PIM kinematics. The results are shown in Figure 7. In the 1PI scheme, there is a sharp 
growth in the average value of the "small" parameter A with increasing /3. The expansion 
parameter in the PIM scheme also increases as a function of (3, but not as quickly. Note 
that this behavior of the partonic expansion parameter A does not translate directly into 
correspondingly large corrections to the threshold expansion for the physical cross sections. 
The singular distributions in s 4 or (1 — z) still enhance the region where A is small, even if the 
integration range covers regions where it is not, and the regions of j3 closer to the endpoint 
are damped by the parton luminosities. However, given the results of the figure, it is not 
surprising that the power corrections are generally larger in 1PI than in PIM kinematics, and 
that they are especially important at the LHC, where the parton luminosities are larger at 
higher values of (5. 

So far, we have focused on the agreement of the corrections in 1PI or PIM kinematics with 
the exact QCD results or with each other. We can gain more information by looking at the 
contributions of the individual terms in the decomposition (39) in the 1PI scheme and also 
at the analogous contributions in the PIM scheme. The assumption of dynamical threshold 
enhancement is that contributions from regions of phase space where the partonic expansion 
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Figure 7: Mean value of the "small" parameter A, where A = (1 — z)/(2y/z) in PIM and 
A = s^/ {2m t ^m1 + s 4 ) in 1PI kinematics. 

parameter A is large, as shown for example in Figure 7, are suppressed due to the properties of 
the PDFs, so one can expect to see a hierarchy between the different terms in the expansion. 
In particular, one would expect that the plus-distributions contribute more than the 5-function 
and of course the power-suppressed contributions contained in R. 

The exact structure of contributions to the total cross section from the different terms are 
shown in Table 2, for the choice fif — rrit- Generally speaking, the P{ distributions are indeed 
enhanced compared to the other terms. In fact, the NNLO contributions from P 3 can be as 
large as the NLO contributions from Pi. However, in both 1PI and PIM kinematics, there 
are large cancellations between the different terms, so that the total NNLO correction turns 
out to be small. In PIM kinematics, this happens at the level of the distributions, and to a 
lesser extent between the 5-function and R terms, which seem to be generically smaller than 
the other terms. In 1PI kinematics, the 5-function terms and especially the power-suppressed 
terms in R are relatively larger than in PIM. We emphasize that the results for the coefficients 
of the Pi distributions are exact, so a full NNLO calculation will change only the 5-function and 
R pieces, moreover in such a way that the cross section from both types of kinematics agrees 
exactly. The numbers above suggest that these terms are relatively small in PIM kinematics, 
because of threshold enhancement, so in order to preserve the good agreement between the 
two types of kinematics they would also need to be small in 1PI kinematics. In that case the 
terms already included in our calculation are the dominant ones at NNLO, although this can 
only be confirmed through the full NNLO results. 

6 Phenomenology 

We now have all the ingredients in place to present detailed phenomenological results for 
the pt and rapidity distributions, the forward-backward asymmetry at the Tevatron, and the 
total cross section. Before doing so, we review some of the necessary inputs. We have already 
described our scheme for evaluating the approximate NNLO corrections in Section 3, and that 



30 









^3 


P2 


Pi 




5 


R 


sum 


Tevatron 


fH 
H 
O 
co 
i — i 

PH 

i — i 




0.57 


0.39 


1.1 

-0.09 


0.30 
-0.19 


-0.06 
-0.37 


-0.38 
-0.31 


0.94 
0.01 


r (1) 


0.39 


0.18 


0.32 
-0.07 


0.08 
-0.14 


-0.11 
-0.12 


-0.06 
-0.12 


0.45 
0.14 


fH 
H 

O 
co 

1 1 

Ph 


r (1) 

r (2) 


0.64 


-0.24 


1.2 

-0.26 


-0.67 
-0.16 


-0.11 
-0.08 


0.20 
0.11 


0.64 
0.02 


r (2) 


0.60 


-0.31 


0.48 
-0.11 


-0.24 
-0.11 


0.12 
-0.08 


0.03 
0.04 


0.38 
0.10 


LHC14 


fH 
H 

O 

CO 

i — i 

Ph 

i — i 


OCT gg' 


8.1 


2.9 


18 
-1.6 


1 

-0.6 


-2 
-5.4 


-9 
-4.6 


7 

-1.1 


r (2) 

fogg 


296 


83 


280 
-20 


14 
-35 


124 
-160 


-125 
-149 


292 
16 


fH 
H 

O 
a> 

1 1 

Ph 


£ ( 2 ) 


5.5 


-2.8 


11 
-1.7 


-9 
-1.0 


-1 
-1.3 


5 

1.1 


6 

-0.1 




287 


-194 


250 
-41 


-189 
-60 


120 
-10 


60 
37 


240 
19 



Table 2: Corrections in pb from the different types of distributions at NLO and NNLO, for 
fi f = m t . 

for evaluating the formulas at threshold in Section 4, where our treatment of PDFs and the 
top-quark mass was also summarized. In addition, we must specify the choice of the matching 
scales jj, s and //^, and of the factorization scale jif. 

We begin by discussing our method for determining an appropriate choice of the soft 
scale fj, s . From theoretical arguments, we expect that the perturbative expansion of the soft 
function should be well-behaved at a scale characteristic of the energy of the real soft radiation, 
which is generally smaller than the hard scales m t and V^- As explained in Section 3.4, a 
general analysis would determine the soft scale by requiring that the corrections from the 
soft function to the double-differential spectrum are well behaved, after integration over the 
partonic variables. We have performed such an analysis for the pt distribution and found 
reasonable results at the Tevatron. This is also the case for the LHC at lower values of pt, 
where the differential cross section is large. As long as we study a relatively modest range in 
Pt, an equally valid procedure for determining the soft scale is to study the corrections to the 
total cross section as a function of /j, s . This automatically samples the regions of phase-space 
where the double-differential cross section is largest. We show the results of such an analysis 
in Figure 8. To isolate the a s correction from the soft function shown there, we pick out the 
piece of the NNLL approximation to the hard scattering kernels arising from evaluate 
the total cross section using only this piece, and divide the result by that at NLL, working 
in the 1PI SCET scheme. We furthermore make the scale choice fj,f — /ih — A* s , which amounts 
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Figure 8: The one-loop correction from the soft function (left) and the hard function (right) 
to the total cross section, divided by the leading-order result. The solid black line is for the 
Tevatron, the dashed red line for LHC7, and the dotted blue line for LHC14. 

to looking at the correction at NLO in fixed-order. As seen from the figure, a well-defined 
minimum in the soft correction appears for fi s ~ 60 GeV at the Tevatron, fi s ~ 80 GeV at the 
LHC with v/i = 7 GeV, and fi s ~ 90 GeV at the LHC with = 14 GeV. We will use these 
as the default choices of fi s in the rest of this section, both for the total cross section and for 
differential distributions. 

One can apply this same procedure to determine an appropriate choice of the hard scale. 
Since the hard function is the same in PIM and 1PI kinematics, we first recall the analysis 
of [55], where it was argued that fih = M is a reasonable default value. Translated to 1PI 
kinematics, this would imply the choice fih = a/§. The actual result for the a s correction 
to the total cross section arising from the hard function as a function of fih is shown in the 
right panel of Figure 8. We isolate this correction as we did for the soft function, except for 
this time we pick out the piece of the NNLL cross section proportional to H^ 1 ', and examine 
the result as a function of fih- As for the analysis with the soft function, we make the choice 
fij = ji h = fi s . At lower values of fih the correction becomes negative and depends strongly on 
the scale. To avoid sensitivity to that region, we will choose fih = 400 GeV by default, which 
is close to the average value of \f§ for the total cross section, and in any case will be varied 
by a factor of two in the error analysis. 

For the factorization scale, we will consider the two different choices fif = m t and fif = 
400 GeV used in [55]. Of course we could also just choose a single, intermediate value of 
fif and vary it in a larger range, but in resummed perturbation theory it is useful to have 
independent variations of the matching scales fih and fi s for the two different values of fif. For 
the differential distributions in the following subsection, on the other hand, we use fif = 2m t 
as the central value. A more refined analysis could use other choices, such as fif = mi for 
the pt distribution, but since we do not study tails of the distributions we prefer to stick to 
a single value which is roughly intermediate between the two values used for the total cross 
section. 
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Figure 9: Dependence of the cross section at NLO+NNLL and approximate NNLO in IPIscet 
on the scales ///, fi s , and /i^. The dot-dashed red lines show the dependence of the cross section 
on fif in approximate NNLO. The dependence of the cross section on the scales at NLO+NNLL 
order is represented by the solid black lines (///), the dashed black lines (/i s ), and the dotted 
black lines (/i/J. 
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We should mention that another method often used to argue for a particular scale choice is 
to look for areas where the scale dependence of the observable is flat. As part of our analysis 
below, we show in Figure 9 the dependence of the total cross section on the scales /i s , fih, and 
/j, f, at NLO+NNLL and approximate NNLO. We note that the scale-dependence of the cross 
section at NLO+NNLL order is indeed flat close to our default values of fi s and /i^, and also 
close to /if — m t . The approximate NNLO results, on the other hand, do not seem to favor a 
particular choice of /if based on this criteria. 

6.1 Rapidity and transverse-momentum distributions 

We now present results for the top-quark rapidity and transverse-momentum distributions. We 
begin by studying rapidity distributions in Figure 10, where we compare results from fixed- 
order and resummed calculations in the 1PI S cet scheme at the Tevatron and the LHC with 
V^s = 7 and 14 TeV. The results within a given perturbative approximation are represented 
as bands indicating the theoretical uncertainties from scale variations. To make the bands in 
fixed-order perturbation theory at a given point in y, we vary the value of the factorization scale 
up and down from its default value at /if — 2m t by a factor of two, and pick out the highest 
and lowest numbers at that point. To make the bands in resummed perturbation theory, we 
vary individually /if, fih, or /i s with the others held fixed at their default values, pick out 
the highest and lowest numbers, and then add the three highest and three lowest numbers in 
quadrature. The results in the figure clearly show that the higher-order corrections contained 
in the approximate NNLO and NLO+NNLL formulas tend to reduce the uncertainties due to 
scale dependence, and slightly raise the central values at a given rapidity to the upper part 
of the fixed-order NLO band. The error bands for the NLO+NNLL and approximate NNLO 
results are of similar size at the Tevatron, but at the LHC the approximate NNLO bands are 
noticeably smaller, a result which we will quantify in more detail when we study the total 
cross section. 

Next, we consider the top-quark transverse- momentum distribution. In this case, we focus 
our analysis on the Tevatron, where experimental measurements are available. In Figure 11 
we show our predictions for this distribution within the different perturbative approxima- 
tions, comparing the fixed-order and resummed results. As with the rapidity distributions, 
the higher-order perturbative corrections serve to decrease the scale dependence, and also to 
slightly raise the central values of the results at a given p T . We compare the NLO+NNLL 
results for the pr distribution with a recent measurement at the Tevatron performed by the 
DO collaboration using the lepton+jets channel [79] in Figure 12, showing also the NLO cal- 
culation for illustration. Since the DO analysis uses m t = 170 GeV, for the purposes of this 
study we deviate from our default choice and also adopt this value. We observe that the slight 
increase in the the differential cross section due to resummation leads to a better agreement 
with the data compared to the NLO predictions. In general, the measured spectrum and the 
NLO+NNLL theory prediction agree within the errors, both in the shape and the normal- 
ization. This is true even at higher values of p?, although one should keep in mind that our 
scale-setting procedure was designed to work in areas where the differential cross section is 
large, and at p T ~ 400 GeV it would also be interesting to perform a more involved analysis 
where ///,/i s , and fj,h are chosen as dynamical functions of px- 
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Figure 10: Left: Fixed-order predictions for the rapidity distribution at LO (light), NLO 
(darker), and approximate NNLO (dark bands) for the Tevatron (top) and LHC (bottom). 
Right: Corresponding predictions at NLL (light) and NLO+NNLL (darker bands) in re- 
summed perturbation theory. The width of the bands reflects the uncertainty of the dis- 
tributions under variations of the matching and factorization scales, as explained in the text. 
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Figure 11: Left: Fixed-order predictions for the pt distribution at LO (light), NLO (darker), 
and approximate NNLO (dark bands) at the Tevatron. Right: Corresponding predictions at 
NLL (light) and NLO+NNLL (darker bands) in resummed perturbation theory. The width 
of the bands reflects the uncertainty of the distributions under variations of the matching and 
factorization scales, as explained in the text. 
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Figure 12: Comparison between our NLO+NNLL predictions, NLO results and recent mea- 
surements from the DO collaboration [79] . The error bands refers to perturbative uncertainties 
related to scale variations. Furthermore we have enlarged the region of bigger pt for better 
comparison. 
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0.2 < Hf/TeV < 0.8 


mt/2 < fif < 2m t 


Aa F B [pb] 


^FB [%] 


Aa FB [pb] 


^FB [%] 


NLL, IPIscet 

NLO leading, 1PI SCET 

NLO 


fl OQ4+0.099 

u.z,o £ ±_ 116 
n 164+0 077 

16S+ 076 
u.±uo_ 048 


4 12+ 2 - 11 

4 40+0-38 
<±.OU_ 32 


0-238™ 
0-2621°;^ 

260+° 140 
u.z,uu_ 084 


4.241»;g 
4-86±g;g 
4.8118^ 


NLO+NNLL, IPIscet 
NNLO approx, IPIscet 


n oqc;+o.o26 
941 +°- 023 


4 61+ - 27 

^• ul -0.26 

4 27+ - 31 

^•^ ' -0.26 


fl ^1 9+0.027 
U.O±^_ 035 

270+ 037 

U.Z, 1 U_ 


4.881g:£ 

4.oi±g : gs 



Table 3: The asymmetric cross section and forward-backward asymmetry at the Tevatron, 
evaluated at different orders in perturbation theory in the pp center-of-mass frame. The errors 
refer to perturbative uncertainties related to scale variations, as explained in the text. 



6.2 The forward-backward asymmetry at the Tevatron 

The top-quark pair charge asymmetry is an important observable at the Tevatron. It originates 
from the difference in the production rates for top and anti-top quarks at fixed scattering angle 
or rapidity [22,23]. Due to charge-conjugation invariance in QCD, it can also be interpreted 
as a forward-backward asymmetry using (9). This asymmetry vanishes at LO in a s and is 
only about 5%-10% at NLO, so it is potentially sensitive to new physics contributions and 
has received much interest for that reason. 

The forward-backward asymmetry depends on the frame in which it is measured. Ex- 
perimentally the asymmetry has been measured both in the laboratory frame and in the tt 
rest frame [7, 8] . In [55] , we used results from PIM kinematics to calculate the asymmetry to 
NLO+NNLL and approximate fixed order in the partonic center-of-mass frame, with the aim 
of improving the NLO fixed-order calculations from [23] and the NLL results from [62]. In 
what follows we use our results for IPI kinematics to obtain NLO+NNLL and approximate 
NNLO predictions directly in the laboratory frame. 

We define the forward-backward asymmetry as the ratio 

A'* B = ^t, (58) 

where the asymmetric cross section for the difference between the production of top quarks in 
the forward and backward directions is 

Acr FB = / dy dp T — — / dy dp T — — - — , (59) 



JO u,y T u,y j _ y+ 



dp T dy J J dp T dy 



with 



U+ 2 l-^l-Am 2 /s T 2 y cosh 2 ?/ s 1 ; 

In Table 3, we show the results for the asymmetric cross section and forward-backward asym- 
metry obtained in different perturbative approximations. In calculating the asymmetry, we 
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first evaluate the numerator and denominator of the ratio A\. B = Aap^/cr to the order indi- 
cated in the table (using our scheme for the PDFs from Table 1), and then further expand 
the ratio itself, as appropriate at that order. We have obtained the fixed-order NLO results 
using the formulas from Appendix A of [23] and cross-checked them using MCFM. As seen 
from the table, the higher-order corrections contained in the NLO+NNLL and approximate 
NNLO formulas stabilize the asymmetric cross section at values not far from the NLO calcu- 
lation with /if = m t . The forward-backward asymmetry shows a wider spread, and its value 
is not correlated with the asymmetric cross section in a straightforward way. For instance, 
at /if — m t , the central value of the asymmetric cross section in fixed order is essentially 
unchanged but has noticeably smaller scale uncertainties at approximate NNLO compared 
to NLO. Yet the asymmetry itself decreases and the errors are not reduced in proportion to 
those for the asymmetric cross section. However, simply taking the higher and lowest numbers 
from the NLO+NNLL and approximate NNLO numbers leads to predictions in the range of 
roughly 4%-5.6%, so the message of this analysis is that the NLO prediction is rather sta- 
ble under radiative corrections. For comparison, the most recent measurement made by the 
CDF collaboration of the inclusive forward-backward asymmetry in the laboratory frame is 
A f FB = (15.0 ± 5.5)% (stat.+sys.) [9]; roughly, the experimental result exceeds the theoretical 
predictions it Table 3 by two standard deviations. The paper [9] also reports a measurement 
of the asymmetry in the tt frame as a function of the pair invariant mass; in the high mass 
region the experimental measurement exceeds the NLO QCD prediction by more than three 
standard deviations. It is possible to calculate the forward-backward asymmetry as a function 
of the invariant mass up to NLO+NNLL accuracy by employing the results of [55]. The results 
of such an analysis will be presented elsewhere [80]. 

6.3 Total Cross Section 

Our main new results concerning the total cross section are the NLO+NNLL and approximate 
NNLO expressions in the IPIscet scheme, obtained here for the first time. In addition to 
presenting these numbers, we also pay attention to a comparison with the PIMscet results, 
and a main outcome of this study is that the total cross sections obtained within the two 
types of kinematics are in good agreement, as long as the sub leading terms identified within 
the SCET formalism are included. This is arguably not the case in the traditional 1PI or PIM 
schemes used in previous calculations, a point we discuss in more detail in Section 6.4. 

Our results for the total cross section within different approximations are summarized 
in Tables 4 and 5. We have also shown the scale dependence of the IPIscet results at 
NLO+NNLL and approximate NNLO in more detail in Figure 9. The main phenomeno- 
logical results are the NLO+NNLL and approximate NNLO numbers in the bottom half of 
the tables, but for comparison we also give the NLO predictions obtained within the differ- 
ent approximations. In that case, we also present the sum of the qq and gg channels alone, 
without the extra piece from the gq and gq channels; this allows for a more direct comparison 
with the leading-singular results in the IPIscet and PIMscet schemes at NLO, neither of 
which includes those pieces. In the fixed-order results, the scale uncertainties are obtained by 
varying the factorization scale up and down by a factor of two. In the resummed results, we 
vary independently fih, fi s and ///up and down by a factor of two with the other scales held 
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Tevatron 


LHC (7TeV) 


LHC (8TeV) 


LHC (14TeV) 


CNLO leading, 1PI S cet 
C"NLO leading, PIM SCET 
0"NLO, qq + gg 
0"NLO 


K QO+0.74+0.33 
>J-»^_0.80-0.22 

c Cfi+0.78+0.31 
U.UU_ 78 _q .20 

5.89±g : n 

c 7Q+0.79+0.33 
°- 1 y -0.80-0.22 


14Q+13+8 
LH:U -16-8 

134±i«±? 
142±# 

1 qq+21+7 

100 -19-7 


91 4+16 + 10 
z *- L4: -22-10 

1 Q9+ 21 + 9 

ls,z, -23-9 

2031- 
1 q9+30+9 


Ocq+35+29 
o,JO -65-30 

7«1 +64+25 
' u± -75-26 

8011?? 
7611^ 5 lf 7 


CNLO+NNLL, IPIscet 
CNLO+NNLL, PIM scet 
0"NNLO approx, IPIscet 
°NNLO approx, PIMscet 


a cq+0.14+0.32 
u - UJ -0. 17-0.23 

a oq+0.19+0.31 
<J.^»_ .20-0.23 

a qn+0.30+0.32 
U.OU_ 3 9_o .23 

a 1 O+0.43+0.31 
u - lz, -0.47-0.23 


157±I X ±1 
149™ 
153™ 
145™ 


99q+9 +io 

zzo -15-ll 

91 9+10+10 
2 iq+2+io 
9(17+11+1° 


8451-1- 

O9fi+40+28 
oz,u -44-29 

847+6+28 
OH: 1 -0-30 

01 1 +38+27 
01 x -25-29 



Table 4: Results for the total cross section in pb, using the default choice /i/ = 400 GeV. The 
first set of errors refers to perturbative uncertainties associated with scale variations, and the 
second to PDF uncertainties. 





Tevatron 


LHC (7TeV) 


LHC (8 TeV) 


LHC (14 TeV) 


CNLO leading, IPIscet 
C"NLO leading, PIMscet 
0"NLO, qq + gg 

0"nlo 


c 7Q+0.20+0.38 
u - ' ^-0. 70-0. 24 

a 49+0.42+0.35 
"•^-O. 76-0. 23 

6.801°;- 

a 79+0.36+0.37 
u - ' ^-0.76-0.24 


1631°^ 
152±I 5 ±1 
160lf 5 

1 cq+20+8 
ida -21-9 


9q9+o +11 

Z " JZ '-14-12 

91 7+ 8 +1° 
^ L 1 -20-11 

2281 6 20 

997+28+11 

^ ' -30-12 


887+o +30 

00 ' -66-32 

8361gl» 

8791- 
889+ 107 + 31 

oo ^-106-32 


CNLO+NNLL, IPIscet 
CNLO+NNLL, PIMscet 
CNNLO approx, IPIscet 
°NNLO approx, PIMscet 


a cc+0.16+0.32 
U.UU_ 14 _ 24 

a 4^+0.18+0.32 
u -^ u -0.19-0.24 

fi fiq+0.00+0.33 
U.UO_ 27-0.24 

a ^9+0.05+0.33 
u - uz, -0.40-0.24 


i5o±;±i 
147+7+8 

LH: ' -6-8 

155™ 
155™ 


91 4+10+10 
^i-^-io-ii 

91 n+10+10 

999+5+11 
z,z,z,_ 3 _ 11 

991 +12+11 

Z,Z,J --12-12 


894+41+28 

oi 1 +45+29 
°- L1 -42-30 

8511fl 2 3? 
860±g±i§ 



Table 5: Same as Table 4, but with the scale choice fij = m t . 



fixed, and then add the three uncertainties together in quadrature to obtain the perturbative 
uncertainties shown in the tables. The resummed results in the PIMscet scheme are obtained 
with the choices of Hh and [x s used in [55]. We have included PDF uncertainties obtained by 
evaluating the cross section with the set of MSTW2008 PDFs at 90% confidence level. 

In all cases, the results in the IPIscet and PIMscet schemes within a given perturbative 
approximation are compatible with one another, once the uncertainties estimated from scale 
variations are taken into account. The IPIscet results tend to come out higher and the 
PIMscet results lower, but the differences are not dramatic. At NLO, the central values of 
the exact results always fall between the predictions obtained by retaining only the leading 
singular pieces in the two types of kinematics. The fif dependence is also taken into account 
relatively well by the two types of kinematics, although at higher collider energies the qg and 
qg channels give an important contribution to the exact NLO result at lower values of /if, 
which is not taken into account by the leading-singular pieces of the threshold expansions. 
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Figure 13: Left: Dependence of the total cross section on top-quark mass defined in the pole- 
mass scheme. The IPIscet and PIMscet NLO+NNLL bands reflect the linearly combined 
scale and PDF uncertainties. The blue band shows the dependence of a DO measurement of 
the total cross section on m t [81]. Right: The same for the LHC but with a comparison to 
recent CMS [11] and ATLAS [12] measurements. 



The NLO+NNLL and approximate NNLO formulas include these NLO corrections through 
the matching, and the fact that these higher-order approximations have a rather small scale 
dependence even without the NNLO corrections from the qg channel can be taken as an 
indication that such corrections are small, but this is a point which should nonetheless be 
kept in mind when applying approximate formulas to the LHC with = 14 TeV. 

Another noticeable pattern has to do not with the difference between the IPIscet and 
PIMscet schemes, but rather between the NLO+NNLL and approximate NNLO results: the 
NLO+NNLL results are higher for /// = 400 GeV, but lower for jj,f = m t . We can learn more 
about the scale variations in the NLO+NNLL and approximate NNLO results in the IPIscet 
scheme by examining Figure 9. Compared to both the NLO+NNLL approximations and the 
approximate NNLO numbers in the PIMscet scheme, the approximate NNLO results in the 
IPIscet scheme have very small scale uncertainties around a given fif, particularly at the 
LHC. When ///is considered in the entire range m t /2 < Hf < 800 GeV, on the other hand, 
the approximate NNLO results are in good agreement with the resummed results in that same 
range. This is not an unreasonable means of comparison, since at a given \x$ the resummed 
results probe scales ranging from /i^ cf ~ 75 GeV to Hh ~ 400 GeV, so in estimating errors in 
fixed-order one should arguably focus on a similar range instead of just using mt/2 < [if < 2mt, 
as is often done in the literature. 

As emphasized in Section 3.4, due to our method for determining the soft scale fi s , the 
NLO+NNLL predictions contain slightly different information than the approximate NNLO 
formulas. In particular, the NNLO expansion of the resummed formulas differ in the structure 
of the P n distributions. For instance, the approximate NNLO results contain P3 distribu- 
tions, but the equivalent terms in the direct NNLO expansion of the resummed formula are 
of the form P2 ln(/z s ////). These contribute at the same order in the counting of RG-improved 
perturbation theory, but there are obviously numerical differences between the two forms of 
the expansion. This issue is discussed in more detail in the Appendix. Since the analysis 
in Figure 8 was done at NLO and at that order one encounters at most Pi distributions, 
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it is worthwhile to ask whether the choice of scale of /j, s deduced there is really appropri- 
ate to account for the mismatch between the two approximations. If this were not the case, 
our NLO+NNLL predictions would become unstable upon the inclusion of the higher-order 
matching corrections and it would make more sense to use the approximate fixed-order NNLO 
results. To address this issue, we have also calculated the cross section where we add the 
NNLO matching coefficients on top of the NLO+NNLL resummation. More precisely, we 
include the pieces of the matching functions specified in (40), but still including the NLL 
evolution matrices in the trace, as in (35). In the case where /if — fj,h — A*s, this approxima- 
tion reduces to the approximate NNLO result (compared to the NLO+NNLL result, which 
reduces to the NLO result in this limit). We have checked that the numerical results in this 
"NNLO+NNLL" approximation are within the uncertainties estimated by the NLO+NNLL 
calculation. This provides evidence that our scale-setting procedure is indeed appropriate for 
effectively including the higher-order corrections. 

In Figure 13 we show the total cross section calculated in the 1PI S cet and PIM S cet schemes 
as a function of the top-quark mass defined in the pole-mass scheme. The bands represent 
linearly combined errors from scale and PDF uncertainties. On the left side we compare these 
with the dependence of a DO measurement on m t at the Tevatron [81]. We have converted 
our error bands to 68% CL (la) to match the confidence level of the experimental errors. On 
the right side we present our predictions for the m r dependent cross section in combination 
with new measurements done by the CMS [11] and ATLAS [12] collaborations. Both analyses 
use a top-quark mass of 172.5 GeV, but we have set them aside in Figure 13 for clarity. 
Both for the Tevatron and LHC there is good agreement between theoretical predictions and 
measurements. It would be interesting to use a more physical mass scheme like MS both in 
theoretical and experimental analyses. This will be discussed in a forthcoming article [80]. 

6.4 Comparison with previous results 

Recently, approximate NNLO cross sections and p T distributions in 1PI kinematics have been 
calculated in [69]. While the resummation formalism leading to these results is different, 
we have checked wherever possible that the structure of the NNLO expansion as written in 
(39) is the same. Therefore, our results in the 1PI scheme (but not the IPIscet scheme, 
which includes a set of corrections in the R(si) term in that equation), should in principle 
be the same. However, a direct comparison is complicated for two main reasons. First, no 
explicit analytic results for the coefficients appearing in (39) were given in [69], which focused 
instead on recalculating the two-loop anomalous dimension matrix obtained in [82-84]. As 
explained in the text after that equation, there are ambiguities in what terms to include in 
the coefficient multiplying the 5-function term. Second, even if the scheme were fixed and the 
analytic expressions agreed, there would still be ambiguities in the numerical implementation, 
as explained in Section 4. In particular, depending on the choice of the variables s', t[ and &[, 
the corrections stemming from the plus-distributions can be different. Since analytic results 
and the choice of scheme were not specified in [69], we cannot compare directly with the 
numerical results given there. 

Instead, we shall present our own results in the 1PI and PIM schemes, and use those to 
compare with the IPIscet and PIM S cet results. For the default scale choice /if — m t — 
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Figure 14: The a s corrections to da/dp for the choice jj,f = m t . The 1PI and IPIscet results 
refer to the scheme used in the present work, while the IPIklmv results uses that from [59]. 
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Tevatron 


LHC (7TeV) 


LHC (14TeV) 


0"NLO leading, 1PI (1PI SC et) 
0"NLO leading, PIM (PIM SC et) 


7.231g;« (6.791°,;?°,) 
6.201°;^ (6.4218:$ 


1831?, (1631°,) 
14311 2 (1521L) 


10241° 7 (8871° 6 ) 
7711° 2 (8361*1) 


°NNLO approx, 1PI (IPIscbt) 
0"NNLO approx, PIM (PIM SC et) 


7.06±g;S (6.631°;°?) 
6.46l°i« (6.621°,;°*) 


18011 (1551^) 
148111 (15511) 


10091^° (8511f) 
8231™ (860lg) 



Table 6: The total cross section in pb in the 1PI and PIM schemes, where the subleading terms 
produced through the SCET expansion are not included. The uncertainties are estimated 
from scale variation in the range mt/2 < [if < 2m t . The results in the IPIscet and PIMscet 
schemes (also shown in Table 5) are listed in parentheses for comparison. 

173.1 GeV, the numbers in these schemes are shown in Table 6. Comparing with Table 5, we 
notice that the numbers in the PIM SCET and PIM schemes are in much better agreement than 
those in the 1PI S cet and 1PI schemes. Actually, the PIM SCET , 1PI SCET , and PIM results are 
compatible with each other within the uncertainties estimated through scale variations, while 
the 1PI results are much higher, especially at the LHC. Since all of our numerical studies 
have indicated that predictions in 1PI kinematics at LHC energies are susceptible to large 
power corrections, we are rather hesitant to give much weight to those particular results. We 
consider the results in PIM kinematics more reliable, and if different kinematics is to be used 
as a criterion for estimating uncertainties, we believe the situation is more accurately reflected 
in Tables 4 and 5 than in Table 6. 

The approximate NNLO numbers in Table 6 are very close to those of [69] for the Tevatron, 
but at the LHC they are roughly 10% higher than those in [69]. While we are unable to pinpoint 
the source of this discrepancy for the reasons explained above, we can compare the leading 
singular terms in our 1PI expansion at NLO with those obtained in [59]. As explained in 
Section 4, these were obtained using a different scheme to evaluate the equivalent analytic 
results at threshold, so the difference between two numbers gives a rough measure of the 
power-suppressed ambiguities inherent to the formulas at threshold. In Figure 14 we compare 
the results of [59] with those obtained in our 1PI scheme, in the 1PI S cet scheme, and with the 
exact results. As seen from the figure, the leading terms in the IPIscet scheme provide the 
best approximation to the full results, and differences due to the choice of scheme get larger as 
the collider energy increases. The numbers in our 1PI scheme are higher than those from [59], 
which may account for the discrepancy with [69]. For reference, the leading singular terms at 
NLO with pif = m t in the scheme of [59] yield cross sections of 7.42 pb at the Tevatron, 174 pb 
at the LHC with yfs = 7 TeV, and 968 pb at the LHC with y/s = 14 TeV. In [69], it was stated 
that in the gg channel the leading singular terms in 1PI kinematics account for over 98% of 
the NLO correction at both the Tevatron and the LHC. For this NLO correction, the exact 
results dX jif = m t are 0.42 pb at the Tevatron and 280 pb at the LHC with y/s = 14 TeV, 
compared to 0.44 pb and 360 pb for the leading singular pieces in the 1PI scheme of [59]. The 
number at the LHC is quite a bit higher than the exact result, and we therefore are not able 
to confirm the statement made in [69], although we cannot rule out that they were obtained 
in yet another (unspecified) scheme for evaluating the formulas at threshold. 
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7 Conclusions 



We have used techniques from SCET to study higher-order perturbative corrections to 1PI 
observables in top-quark pair production at hadron colliders. In particular, we have calcu- 
lated and collected the perturbative ingredients needed for an analysis at NLO+NNLL and 
approximate NNLO, and applied them to transverse- momentum and rapidity distributions 
as well as the total production cross sections at the Tevatron and LHC, and to the forward- 
backward asymmetry at the Tevatron. The details of factorization and resummation within 
our approach were mainly taken over from the analogous study in PIM kinematics presented 
in [55]. However, differences arise from the phase-space integrations for real emission in the 
soft limit. They affect the matrix- valued soft matching function, which must be calculated 
from scratch. We presented results for this matching function at NLO in fixed-order pertur- 
bation theory, which were previously unknown, and also addressed the differences in its RG 
equation needed for the resummation of soft gluon corrections to all orders. Along with the 
anomalous dimensions and matrix- valued hard matching functions already known from [55], 
these are the essential ingredients for an NLO+NNLL calculation. 

Even apart from the new perturbative information included in our results, an important 
aspect of the SCET approach to threshold resummation is its ability to include a class of loga- 
rithmic corrections, which are formally subleading in the exact threshold limit but nonetheless 
give non-negligible numerical contributions. These arise in the SCET formalism by noticing 
that in a suitable reference frame the single dimensionful parameter characterizing soft radia- 
tion is the total soft energy E s of the emitted gluons, and expressing the logarithmic corrections 
in terms of Yn.(E s /n), without further expansion. The subleading terms generated in this way 
are of the form ln(l + Si/mf)/si in 1PI kinematics, and lnz/(l — z) in PIM kinematics (with 
z = M 2 /s). We have implemented these new terms in our numerical analysis summarized 
below, where results including them are referred to as the 1PI SCET and PIM SCET schemes. 
Results neglecting them are referred to as the 1PI and PIM schemes and should correspond 
with previous results at a given order in perturbation theory. 

Our improved calculations of the hard-scattering kernels near threshold are a priori only 
important in regions of phase space where kinematic restrictions imply the partonic threshold 
limit S4 — > 0. However, the differential cross section in such regions is very small. Beyond 
that, the improved predictions are important if a dynamical threshold enhancement occurs, 
due to a sharp fall-off of the PDFs away from integration regions where s 4 — > 0. In Section 5 
we performed a very detailed discussion of the mechanism of threshold enhancement, and 
compared its effectiveness in 1PI and PIM kinematics at the level of differential distributions 
and the total cross section. The main finding was that the leading power corrections to the 
threshold limit in 1PI kinematics, which scale as E s /m t ~ s^/mf, are in general larger than 
the leading corrections in PIM kinematics, which scale as E s /M ~ (1 — z). At the LHC with 
^J~s = 14 TeV these corrections can become quite large. This makes us somewhat reserved 
about the utility of the threshold expansion in 1PI kinematics, especially at larger values of pt- 
Including the extra subleading terms unique to the IPIscet approximation seemed to account 
for the power corrections to a good degree of accuracy in the qq channel and also to some 
level in the gg channel, at least in regions of phase space where the differential cross section 
is large. Since the top-quark pair production cross section at the Tevatron is dominated by 
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the quark-antiquark channel, the predictions obtained in the IPIscet scheme at the Tevatron 
center-of-mass energy are expected to provide a reliable approximation to the exact results. 
In all cases, but especially at the LHC with y/s = 14 TeV, the pure 1PI expansion gave larger 
corrections compared to any other approximation, and we are rather skeptical of any results 
obtained in that scheme. For the reasons mentioned above, we limited the phenomenological 
analysis to areas where the differential distributions are fairly large, or else to integrated 
quantities such as the forward-backward asymmetry and the total cross section. Our results 
for the rapidity and transverse- momentum distributions can be found in Figures 10 and 11. 
As expected, including the higher-order perturbative corrections in the form of NLO+NNLL 
or approximate NNLO predictions leads to a stabilization of the results under variations of 
the different perturbative scales. Our main results for the total cross section can be found in 
Tables 4 and 5. Although we consider the PIMscet numbers more trustworthy, it is reassuring 
that the IPIscet numbers are actually in very good agreement with them. This is not the 
case for the pure 1PI numbers, which come out significantly larger in all cases. If different 
types of kinematics is to be used as a source of uncertainty in predictions based on soft-gluon 
resummation, we believe the true results are more accurately reflected in Tables 4 and 5. 

In Section 6.3 and the Appendix we analyzed in detail the origin of the numerical dif- 
ferences between NLO+NNLL and approximate NNLO calculations. The structure of the 
plus-distributions considered in the two approaches is different, since in the resummation ap- 
proach plus-distributions and logarithms of scale ratios are formally of the same order. We 
showed why our choice of the soft and hard scales for the NLO+NNLL calculations is sound 
and how the resummed results are stable with respect to reasonable variations of these scales. 

In conclusion, we obtained NLO+NNLL predictions for the top-pair transverse-momentum 
and rapidity distributions at hadron colliders. In the case of the transverse-momentum distri- 
bution at the Tevatron, the predictions presented here can be already compared with experi- 
mental measurements. By integrating the results for the differential distributions we obtained 
predictions for fully inclusive observables such as the forward-backward asymmetry and the 
total cross section. For the total cross section in particular, we analyzed the relations among 
the IPIscet results obtained here, the PIM S cet results of [55], and the ones obtained with a 
different resummation approach [69]. All the evidence gathered shows that the SCET-based 
approach provides consistent and reliable results in both types of kinematics. 
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A Appendix 



In our numerical analysis we considered two different ways of including higher-order perturba- 
tive corrections to the hard-scattering kernels: constructing approximate expansions to NNLO 
in fixed order, and all-orders resummmation to NNLL accuracy. We emphasized that for a 
fixed numerical value of fi s , the truncation of the NNLL series to NNLO in a s contains a 
different structure of corrections than the approximate NNLO formula. Here we explain this 
statement in more detail, and show explicitly what types of higher-order corrections the mas- 
ter formula (35) resums. To do so, it is sufficient to ignore the matrix structure of the RG 
equations, and just consider hard and soft functions which are simple functions of their argu- 
ments. Then the treatment of the hard function is straightforward, and the complication for 
the soft function is its non-locality. 

To explain the issue, it is simplest to start with the hard function. Ignoring its matrix 
structure and dependence on kinematic invariants, its RG equation is of the form 



/i^H(L h ,a s (ii)) = ( 2r cusp 



hi— + 2>A //( /.,, n J/<)). 



(61) 



where Lh = ln(M 2 //i 2 ). The parameter M is to be understood as a generic hard scale; in 
the specific case of 1PI kinematics, it would be M = V§. The RG equation can be used to 
generate higher-order terms in the perturbative expansion of H . For instance, we can solve 
the equation as a fixed-order series in a s using the ansatz 



H(L h ,a s ((j,)) = a 2 s (fi) 
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n=0 



(62) 



In terms of the lowest-order coefficient and the anomalous dimensions and ^-function, whose 
expansions we define as 



cusp 



a 



i — i ( — i • (^V • ... . 



47r \4irJ ' ~ " \AtiJ 



/3(a s ) = -2a s 



n a s n / a s \ 2 „ /o s \ 3 



(63) 



and similarly for 7 h , the expansion coefficients at NLO can be constructed as 



H(L h ,a s (fi)) = a 2 s (fi)\h^ + 



An 



+ . . 



(64) 

One can obviously generalize this to any order in a s and calculate the coefficients of the 
logarithms at a given order in terms of the anomalous dimensions and lower-order matching 
coefficients. This is the method used in constructing approximate fixed-order expansions. 

In the effective-theory analysis, one assumes the presence of a second, widely separated 
scale fi s M and uses the counting of RG-improved perturbation theory, i.e. ln(/i s /M) ~ 
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l/a s . Then the higher-order corrections contain large logarithms, which can be resummed to 
all-orders by using the exact solution to the RG equation. This solution reads [85] 

(M 2 \ / M\~ 2rih 

\n-t-,aM) [—) , (65) 

where r/h = 2ar(/i/i, /x). The Sudakov exponent and normal anomalous exponent are 

S^) = - I da^^l J a r (^) = - J da , (66) 

and similarly for a 7 h. 

The all-orders solution does not actually depend on /ih, as indicated by the notation. The 
same is true if the matching coefficient and the exponentials are consistently re-expanded as 
a series in a s (/x) in fixed order, in which case one just gets back the approximate formulas 
above. In practice, however, one must truncate the result at a given level of accuracy (e.g. 
NNLL), and beyond that level a residual dependence on fih remains. To avoid large logarithms 
in the matching coefficients, one chooses jih ~ M and runs to the scale /if using the all-orders 
solution. Then the exponential factors resum logarithms which count as \n(M 2 /fx 2 ) and are 
large for /if ~ fj, s . To see this explicitly, we get rid of a s (/ih) everywhere by using (with 
a(fi f ) = a s (fi f )/4n) 



aW = ^- %^lnX + ...; X = l + A) a M In 4 , (67) 



and re-expand the solution in (65) as a series in a(/j,f). For reference, the expansion of the 
Sudakov factor and anomalous exponent to NNLO read 



(68) 



S(n h ,n f ) = -a(/i/)y L 2 hf + a(fi f ) 2 (^^L 3 hf -^L) 

ar(Hh, Vf) = a(fif) y L hf + a(fi f ) 2 (^-^-L 2 hf + + . . . , (69) 

where L h f = In^l/ fj, 2 ). The expansion of a 7 is identical to the one in (69), with the replace- 
ments Ti — > 7j. 

Using the above equations to expand (65) to NLO, one recovers the NLO solution (64). 
All the dependence on fih drops out to that order, as long as we keep the one-loop matching 
correction. However, this is not the case if we expand our NLL approximation of the resummed 
hard function to NLO. In that case, the NLO matching coefficient is of higher-order in the 
counting and not included in the formula, so after expansion to NLO dependence on the scale 
\ih remains. The direct expansion of our NLL formula at NLO in fixed order reads 



a 2 s ^)a^ f ){-h^\A 



^ln^ + r ln^- + ( 7 * + 2/9b) 



(70) 
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This is necessarily different than the "approximate NLO" formula one would deduce by drop- 
ping the coefficient h^ 1,0 ^ from (64), because it depends on fj,h, but if we set = M it is the 
same, and it is for this reason that one can still say the NLL solution "resums logarithms of the 
form \n(M//if) to all orders", although a more accurate statement would be that it "resums 
logarithms of the form ln(/x/ l //xj) to all orders", which includes the possibility of other choices 
such as fih — m t- Given this fact, it makes little sense to construct an approximate formula 
for a quantity such as the hard function: if there are large logarithms, it is just as easy to sum 
them to all orders as it is to construct the fixed-order expansion, and if the logarithms are not 
large, there is no reason to include that subset of the higher-order corrections without the full 
answer. 

We can repeat the analysis above to compare the structure of approximate fixed-order 
expansions and resummed formulas for the soft function. In this case the RG equation is non- 
local, and to solve for the momentum-space soft function one uses the technique of Laplace 
transforms [71]. The solution for the resummed momentum-space soft function is 



S(co,fi f ) = e -^"^ + ^^'^s{d rn n a ) - — 

to \(J, S 



u \ 2r> e~ 2 ^ 



r(2n) ' 



(71) 



where rj = 2ar(/x s , fJ-f), and s is the Laplace-transformed function, which satisfies the local RG 
equation 



2T, 



M 2 



CUS P ..2 
fj, 



M 2 



In — + 27 s s In— ,a s (/i) . 



(72) 



In this case, approximate formulas in fixed order are obtained by first constructing the solution 
to s using the local RG equation. To NNLO, we use the ansatz 



~/ T „ f,,\\ i | a s(fj) V^„(l,n)rn i f a s(^) \ ST^ J2,n) rn , 

s(L,a s (fi)) = 1 + -^- 2^ s L + l^r) L +•••' 

n=0 ^ ' n=0 

where we set s*- ' ^ = 1 for simplicity. The explicit solution to NNLO reads 



(73) 



s(L,a s (n)) = 1 + 
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+ 



2 J 



+ (7i-/3os (1 ' 0) +7oS (1 ' 0) )^ + s (2 ' 0) 



(74) 



To turn this into an approximate NNLO formula for the momentum-space soft function 
S(u>,/if), one must take the limit /i s — /if and derive replacement rules analogous to (38). 
This is readily done using the expansion 



-(-) =±5(u) + J2-D n (u) V \ 



(75) 
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where the D„ are defined as 



D n {u) 



ln n 



U! 



(If 



(76) 



As was the case with the hard function, the all-orders solution for the resummed soft 
function does not actually depend on the scale /i s , but its truncation to a given logarithmic 
order (e.g. NNLL) introduces residual scale dependence. As explained earlier in the paper, 
our method is to choose /x s to be close to the numerical value where the corrections from the 
soft function to the (differential) cross section are minimal. We then adopt the parametric 
counting ji s ~ oj and apply RG-improved perturbation theory with In /! s //i/ ~ l/ai s , and 
the exponential factors resum logarithms of the form In fi s /fif to all orders. Since the scale 
fi s is dynamically generated through the numerical analysis, it does not appear in the fixed- 
order calculation, so the resummation formula deals with different types of corrections than 
the approximate fixed-order formulas. 4 For this reason, the structure of D n distributions 
appearing at a given order is not the same in the two approaches. 

We now show this in more detail, working first to NLO. In this case, the "approximate 
NLO" formula derived from the solution (74) reads 



S(ijj,Hf) « 1 + a (/•«/) 
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AToD^uj) + 2Y D (u) - —T 5(cu) 



(77) 



This should be compared with the expansion of the NLL formula to NLO in fixed order, for 
which the NLO correction reads (with L s = In /J%/n 2 f) 



a(fi f ) 



2J J- Do' . (^IV.; • j s L s )5(uj) 



+ . . . . 



(78) 



The "approximate NLO" formula has D\ distributions, while the NLO expansion of our NLL 
formula has only D distributions. In the counting of RG-improved perturbation theory, 
however, D ln(/x s //Xf) ~ Di, so the tower of logarithms produced by the expansion of the 
NLL formulas is of course correct. The analogous formula at NNLO is rather lengthy, but to 
illustrate its structure, we focus on the leading correction in the logarithmic power counting, 
which reads 

12L s D 2 (uj) - QLlD^u) + L 3 S D - ^L 4 s 5(u) , (79) 



a{/i f ) 



2 r 2 
1 



while the leading term of the "approximate NNLO" formula is 

a(pif) 2 SrlD 3 (uj). 



(80) 



Again, for oj ~ /i s the terms in the two equations are of the same parameteric order but contain 
different types of distributions: the resummed formulas generate at most D 2 distributions, 
while the approximate NNLO formulas generate D 3 distributions. 

From the discussion above it should be obvious that our formula does not literally resum the 
highest tower of D n distributions to all orders, but rather terms which count that way in RG- 
improved perturbation theory, when the dynamically generated soft scale satisfies fi s ~ u. In 



4 In [65] it was shown that the soft scale decreases as the PDFs fall off more quickly away from values of x 
where uj ~ 0, so the formulas effectively resum logarithms of the slopes of the PDFs. 
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the case of the hard function, we noted that the fixed-order corrections produced by expanding 
the resummed formula were equal to those in the approximate formula for the special choice 
Hh = M. For the soft function, there is no numerical value of ji s for which this would be true, 
but the two are equal if we replace L s — » d v and take the derivatives before re-expanding r\ in 
a s (/if). This procedure can be generalized to all orders by evaluating the formula 




2r ? e ~2lEV 



T(2 V ) 



(81) 



where the factor in the curly brackets is understood to be expanded to all orders as a series in 
a s (/j,f) and so is a function only of ln(/i^/ /i|). In this way, we exponentiate the derivatives with 
respect to rj, which are what generate the highest-order distributions, and the expansion of 
the above formula to any given accuracy in fixed-order reproduces the approximate formulas. 
For instance, if we include the exact one-loop matching (which is just s^ 1,0 )) and the two-loop 
anomalous dimensions, the expansion of the above object gives back our approximate NNLO 
formula for the soft function, plus higher-order terms that resum all the higher-order D n dis- 
tributions at NNLL order, after converting the derivatives with respect to rj with replacement 
rules. While this procedure generalizes choosing fj,^ = M in the hard function, it is by no 
means the same conceptually. The exact hard function is independent of fit, so varying it 
around values fih ~ M gives a way of estimating the higher-order terms. However, it would 
make no sense to replace, for instance, ln(/Xg//x^) — > Cq8 v in (81), with Cq 7^ 1, since the 
derivatives generate both /x-independent and /x-dependent terms. 

The conclusion of this discussion is that results based on approximate NNLO formulas 
contain different information than those based on NLO+NNLL resummation. This is not just 
due to a truncation of the NNLL series to NNLO, but also to the fact that the resummation 
formula exponentiates logarithms depending on the ratio In fi s /Hf in combination with higher- 
order logarithmic plus-distributions, with /x s a dynamically generated numerical soft scale. 
Such logarithms do not appear in the fixed-order calculation, which is independent of /i s . 
Therefore, the choice between using approximate NNLO and NLO+NNLL amounts to whether 
one takes seriously the improved convergence of the soft function at a numerically small soft 
scale. If so, one should use the resummed formulas, if not, one should use the approximate 
fixed-order calculations. In practice, this question can only be answered after a numerical 
analysis, and for this reason, we have evaluated both types of formulas in the studies in 
Sections 5 and 6. 



50 



References 

[1] CDF and DO Collaboration, arXiv:1007.3178 [hep-ex]. 

[2] H. Flacher, M. Goebel, J. Haller, A. Hocker, K. Monig and J. Stelzer, Eur. Phys. J. C 
60, 543 (2009) [arXiv:0811.0009 [hep-ph]]. 

[3] R. Frederix and F. Maltoni, JHEP 0901, 047 (2009) [arXiv:0712.2355 [hep-ph]]. 

[4] V. Barger, T. Han and D. G. E. Walker, Phys. Rev. Lett. 100, 031801 (2008) [arXiv:hep- 
ph/0612016]. 

[5] U. Baur and L. H. Orr, Phys. Rev. D 77, 114001 (2008) [arXiv:0803.1160 [hep-ph]]. 

[6] C. T. Hill and S. J. Parke, Phys. Rev. D 49, 4454 (1994) [arXiv:hep-ph/93 12324]. 

[7] T. Aaltonen et al. [CDF Collaboration], Phys. Rev. Lett. 101, 202001 (2008) 
[arXiv:0806.2472 [hep-ex]]. 

[8] V. M. Abazov et al. [DO Collaboration], Phys. Rev. Lett. 100, 142002 (2008) 
[arXiv:0712.0851 [hep-ex]]. 

[9] T. Aaltonen et al. [CDF Collaboration], arXiv: 1101. 0034 [hep-ex]. 

[10] F. Abe et al. [CDF Collaboration], Phys. Rev. Lett. 74, 2626 (1995) [arXiv:hep- 
ex/9503002]. 

[11] V. Khachatryan et al. [CMS Collaboration], Phys. Lett. B 695, 424 (2011) 
[arXiv: 1010.5994 [hep-ex]]. 

[12] G. Aad et al. [ATLAS Collaboration], arXiv: 1012.1792 [hep-ex]. 

[13] W. Bernreuther, arXiv:1008.3819 [hep-ph]. 

[14] W. Bernreuther, J. Phys. G 35, 083001 (2008) [arXiv:0805.1333 [hep-ph]]. 

[15] P. Nason, S. Dawson and R. K. Ellis, Nucl. Phys. B 303, 607 (1988). 

[16] W. Beenakker, H. Kuijf, W. L. van Neerven and J. Smith, Phys. Rev. D 40, 54 (1989). 

[17] W. Beenakker, W. L. van Neerven, R. Meng, G. A. Schuler and J. Smith, Nucl. Phys. B 
351, 507 (1991). 

[18] M. Czakon and A. Mitov, Nucl. Phys. B 824, 111 (2010) [arXiv:0811.4119 [hep-ph]]. 

[19] P. Nason, S. Dawson and R. K. Ellis, Nucl. Phys. B 327, 49 (1989) [Erratum-ibid. B 335, 
260 (1990)]. 

[20] M. L. Mangano, P. Nason and G. Ridolfi, Nucl. Phys. B 373, 295 (1992). 



51 



[21] S. Frixione, M. L. Mangano, P. Nason and G. Ridolfi, Phys. Lett. B 351, 555 (1995) 
[arXiv:hep-ph/9503213]. 

[22] J. H. Kiihn and G. Rodrigo, Phys. Rev. Lett. 81, 49 (1998) [arXiv:hep-ph/9802268]. 

[23] J. H. Kiihn and G. Rodrigo, Phys. Rev. D 59, 054017 (1999) [arXiv:hep-ph/9807420]. 

[24] M. Czakon, A. Mitov and S. Moch, Phys. Lett. B 651, 147 (2007) [arXiv: 0705. 1975 
[hep-ph]]. 

[25] M. Czakon, A. Mitov and S. Moch, Nucl. Phys. B 798, 210 (2008) [arXiv:0707.4139 
[hep-ph]]. 

[26] M. Czakon, Phys. Lett. B 664, 307 (2008) [arXiv:0803.1400 [hep-ph]]. 

[27] R. Bonciani, A. Ferroglia, T. Gehrmann, D. Maitre and C. Studerus, JHEP 0807, 129 
(2008) [arXiv:0806.2301 [hep-ph]]. 

[28] R. Bonciani, A. Ferroglia, T. Gehrmann and C. Studerus, JHEP 0908, 067 (2009) 
[arXiv:0906.3671 [hep-ph]]. 

[29] R. Bonciani, A. Ferroglia, T. Gehrmann, A. Manteuffel and C. Studerus, JHEP 1101, 
102 (2011) [arXiv:1011.6661 [hep-ph]]. 

[30] J. G. Korner, Z. Merebashvili and M. Rogal, Phys. Rev. D 77, 094011 (2008) 
[arXiv:0802.0106 [hep-ph]]. 

[31] C. Anastasiou and S. M. Aybat, Phys. Rev. D 78, 114006 (2008) [arXiv:0809.1355 [hep- 
ph]]. 

[32] B. Kniehl, Z. Merebashvili, J. G. Korner and M. Rogal, Phys. Rev. D 78, 094013 (2008) 
[arXiv:0809.3980 [hep-ph]]. 

[33] S. Dittmaier, P. Uwer and S. Weinzierl, Phys. Rev. Lett. 98, 262002 (2007) [arXiv:hep- 
ph/0703120]. 

[34] S. Dittmaier, P. Uwer and S. Weinzierl, Nucl. Phys. Proc. Suppl. 183, 196 (2008) 
[arXiv:0807.1223 [hep-ph]]. 

[35] G. Bevilacqua, M. Czakon, C. G. Papadopoulos and M. Worek, Phys. Rev. Lett. 104, 
162002 (2010) [arXiv:1002.4009 [hep-ph]]. 

[36] K. Melnikov and M. Schulze, Nucl. Phys. B 840, 129 (2010) [arXiv: 1004.3284 [hep-ph]]. 

[37] M. Czakon, Phys. Lett. B 693, 259 (2010) [arXiv: 1005.0274 [hep-ph]]. 

[38] M. Czakon, arXiv: 1101. 0642 [hep-ph]. 

[39] G. Sterman, Nucl. Phys. B 281, 310 (1987). 



52 



[40] S. Catani and L. Trentadue, Nucl. Phys. B 327, 323 (1989). 

[41] E. Laenen, J. Smith and W. L. van Neerven, Nucl. Phys. B 369, 543 (1992). 

[42] E. Laenen, J. Smith and W. L. van Neerven, Phys. Lett. B 321, 254 (1994) [arXiv:hep- 
ph/9310233]. 

[43] E. L. Berger and H. Contopanagos, Phys. Lett. B 361, 115 (1995) [arXiv:hep-ph/9507363]. 

[44] E. L. Berger and H. Contopanagos, Phys. Rev. D 54, 3085 (1996) [arXiv:hep-ph/9603326]. 

[45] E. L. Berger and H. Contopanagos, Phys. Rev. D 57, 253 (1998) [arXiv:hep-ph/9706206]. 

[46] S. Catani, M. L. Mangano, P. Nason and L. Trentadue, Phys. Lett. B 378, 329 (1996) 
[arXiv:hep-ph/9602208] . 

[47] R. Bonciani, S. Catani, M. L. Mangano and P. Nason, Nucl. Phys. B 529, 424 (1998) 
[Erratum-ibid. B 803, 234 (2008)] [arXiv:hep-ph/9801375]. 

[48] S. Moch and P. Uwer, Phys. Rev. D 78, 034003 (2008) [arXiv:0804.1476 [hep-ph]]. 

[49] M. Czakon and A. Mitov, Phys. Lett. B 680, 154 (2009) [arXiv:0812.0353 [hep-ph]]. 

[50] U. Langenfeld, S. Moch and P. Uwer, Phys. Rev. D 80, 054009 (2009) [arXiv:0906.5273 
[hep-ph]]. 

[51] M. Beneke, P. Falgari and C. Schwinn, Nucl. Phys. B 828, 69 (2010) [arXiv:0907.1443 
[hep-ph]]. 

[52] M. Czakon, A. Mitov and G. Sterman, Phys. Rev. D 80, 074017 (2009) [arXiv:0907.1790 
[hep-ph]]. 

[53] M. Beneke, M. Czakon, P. Falgari, A. Mitov and C. Schwinn, Phys. Lett. B 690, 483 
(2010) [arXiv:0911.5166 [hep-ph]]. 

[54] M. Beneke, P. Falgari and C. Schwinn, Nucl. Phys. B 842, 414 (2011) [arXiv:1007.5414 
[hep-ph]]. 

[55] V. Ahrens, A. Ferroglia, M. Neubert, B. D. Pecjak and L. L. Yang, JHEP 1009, 097 
(2010) [arXiv: 1003.5827 [hep-ph]]. 

[56] N. Kidonakis and G. Sterman, Phys. Lett. B 387, 867 (1996). 

[57] N. Kidonakis and G. Sterman, Nucl. Phys. B 505, 321 (1997) [arXiv:hep-ph/9705234]. 

[58] N. Kidonakis, Phys. Rev. D 64, 014009 (2001) [arXiv:hep-ph/0010002]. 

[59] N. Kidonakis, E. Laenen, S. Moch and R. Vogt, Phys. Rev. D 64, 114001 (2001) 
[arXiv:hep-ph/0105041]. 



53 



[60] N. Kidonakis and R. Vogt, Phys. Rev. D 68, 114014 (2003) [arXiv:hep-ph/0308222]. 

[61] A. Banfi and E. Laenen, Phys. Rev. D 71, 034003 (2005) [arXiv:hep-ph/0411241]. 

[62] L. G. Almeida, G. Sterman and W. Vogelsang, Phys. Rev. D 78, 014008 (2008) 
[arXiv: 0805. 1885 [hep-ph]]. 

[63] N. Kidonakis and R. Vogt, Phys. Rev. D 78, 074005 (2008) [arXiv:0805.3844 [hep-ph]]. 

[64] V. Ahrens, A. Ferroglia, M. Neubert, B. D. Pecjak and L. L. Yang, Phys. Lett. B 687, 
331 (2010) [arXiv:0912.3375 [hep-ph]]. 

[65] T. Becher, M. Neubert and G. Xu, JHEP 0807, 030 (2008) [arXiv:0710.0680 [hep-ph]]. 

[66] V. Ahrens, T. Becher, M. Neubert and L. L. Yang, Phys. Rev. D 79, 033013 (2009) 
[arXiv:0808.3008 [hep-ph]]. 

[67] V. Ahrens, T. Becher, M. Neubert and L. L. Yang, Eur. Phys. J. C 62, 333 (2009) 
[arXiv: 0809. 4283 [hep-ph]]. 

[68] E. Laenen, G. Oderda and G. Sterman, Phys. Lett. B 438, 173 (1998) [arXiv:hep- 
ph/9806467]. 

[69] N. Kidonakis, Phys. Rev. D 82, 114030 (2010) [arXiv:1009.4935 [hep-ph]]. 

[70] S. Fleming, A. H. Hoang, S. Mantry and I. W. Stewart, Phys. Rev. D 77, 074010 (2008) 
[arXiv:hep-ph/0703207]. 

[71] T. Becher and M. Neubert, Phys. Rev. Lett. 97, 082001 (2006) [arXiv:hep-ph/0605050]. 

[72] T. Becher, M. Neubert and B. D. Pecjak, JHEP 0701, 076 (2007) [arXiv:hep-ph/0607228]. 

[73] C. W. Bauer, N. D. Dunn and A. Hornig, arXiv: 1010.0243 [hep-ph]. 

[74] A. D. Martin, W. J. Stirling, R. S. Thorne and G. Watt, Eur. Phys. J. C 64, 653 (2009) 
[arXiv:0905.3531 [hep-ph]]. 

[75] H. L. Lai, J. Huston, Z. Li, P. Nadolsky, J. Pumplin, D. Stump and C. P. Yuan, Phys. 
Rev. D 82, 054021 (2010) [arXiv: 1004.4624 [hep-ph]]. 

[76] J. M. Campbell and R. K. Ellis, Phys. Rev. D 62, 114012 (2000) [arXiv:hep-ph/0006304]. 

[77] J. Alwall et al, JHEP 0709, 028 (2007) [arXiv:0706.2334 [hep-ph]]. 

[78] R. Frederix, S. Frixione, F. Maltoni and T. Stelzer, JHEP 0910, 003 (2009) 
[arXiv: 0908. 4272 [hep-ph]]. 

[79] V. M. Abazov et al. [DO Collaboration], Phys. Lett. B 693, 515 (2010) [arXiv: 1001. 1900 
[hep-ex]]. 



54 



[80] V. Ahrens, A. Ferroglia, M. Neubert, B. D. Pecjak and L. L. Yang, in preparation. 

[81] V. M. Abazov et al. [DO Collaboration], Phys. Rev. D 80, 071102 (2009) [arXiv:0903.5525 
[hep-ex]]. 

[82] T. Becher and M. Neubert, Phys. Rev. D 79, 125004 (2009) [Erratum-ibid. D 80, 109901 
(2009)] [arXiv:0904.1021 [hep-ph]]. 

[83] A. Ferroglia, M. Neubert, B. D. Pecjak and L. L. Yang, Phys. Rev. Lett. 103, 201601 
(2009) [arXiv:0907.4791 [hep-ph]]. 

[84] A. Ferroglia, M. Neubert, B. D. Pecjak and L. L. Yang, JHEP 0911, 062 (2009) 
[arXiv:0908.3676 [hep-ph]]. 

[85] S. W. Bosch, R. J. Hill, B. O. Lange and M. Neubert, Phys. Rev. D 67, 094014 (2003) 
[arXiv:hep-ph/0301123]. 



55 



